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Abstract 

A summary of the constraints from the ATLAS experiment on 7?-parity-conserving super- 
symmetry is presented. Results from 22 separate ATLAS searches are considered, each 
based on analysis of up to 20.3 fb”' of proton-proton collision data at centre-of-mass ener¬ 
gies of ^/s = 7 and 8 TeV at the Large Hadron Collider. The results are interpreted in the 
context of the 19-parameter phenomenological minimal supersymmetric standard model, in 
which the lightest supersymmetric particle is a neutralino, taking into account constraints 
from previous precision electroweak and flavour measurements as well as from dark matter 
related measurements. The results are presented in terms of constraints on supersymmetric 
particle masses and are compared to limits from simplified models. The impact of ATLAS 
searches on parameters such as the dark matter relic density, the couplings of the observed 
Higgs boson, and the degree of electroweak fine-tuning is also shown. Spectra for surviving 
supersymmetry model points with low fine-tunings are presented. 
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1. Introduction 

During the first run of the Large Hadron Collider (LHC), the ATLAS [1] and CMS [2] Collaborations 
performed a wide range of searches for supersymmetry (SUSY) [3-11], using proton-proton {pp) colli¬ 
sion data at centre-of-mass energies ( ^/s) of 7 and 8 TeV. SUSY, a theoretically favoured framework for 
extending the Standard Model (SM), is able to address some of its unanswered questions, particularly the 
hierarchy problem [12-15], which is related to the fine-tuning needed to obtain the correct mass for the 
observed Higgs boson. SUSY can also provide credible dark matter candidates [16, 17] and can improve 
the unification of the electroweak and strong interactions [18-26]. 

The minimal supersymmetric extension of the Standard Model (MSSM) [27-31] predicts partners for 
each of the SM states. It predicts a pair of scalar partners - one for each fermion chirality - for each 
of the SM quarks and leptons. These spin-zero partner particles are known as squarks (q) and sleptons 
{{) respectively. In the first two generations the pair of chiral partners is largely unmixed, so the mass 
states can be labelled and Sr, where the L and R subscripts denote the scalar partners of the left- and 
right-handed Standard Model fermion states respectively. In the third generation of quarks and leptons 
the mixing between the scalars is larger, and the mixed states are labelled by their mass indices e.g. t\ 
and t 2 , where is lighter by construction. Each state in the SM gluon colour octet has a spin-half partner 
known as a gluino g. There are a total of eight spin-half partners of the electroweak gauge and Higgs 
bosons: the neutral bino (superpartner of the U(l) gauge field); the winos, which are a charged pair and a 
neutral particle (superpartners of the W bosons of the SU(2 )l gauge fields); and the Higgsinos, which are 
two neutral particles and a charged pair (superpartners of the Higgs field’s degrees of freedom). The bino, 
winos and Higgsinos mix to form four charged states called charginos^* 2 > four neutral states known 
as neutralinos (where the index i lies in the range 1 to 4, ordered by increasing neutralino mass). The 
charginos and neutralinos are collectively referred to as electroweakinos. 

Since no statistically significant signals consistent with supersymmetry have yet been observed at the 
LHC, searches have been used to constrain the allowed supersymmetric model space. In the case of 
searches for supersymmetry, this typically results in setting lower limits on the masses of the pertinent 
supersymmetric partner particles (sparticles). 

This paper presents the combined sensitivity and constraints from 22 separate ATLAS analyses of the 
Run 1 LHC dataset, using centre-of-mass energies of 7 and 8 TeV and an integrated luminosity of up to 
20.3 fb“^ Direct searches for the decay products of the sparticles listed above are considered together 
with searches for disappearing tracks, long-lived charged particles, monojet signatures, and a dedicated 
search for the heavier neutral Higgs bosons, also expected in the MSSM. More details about the searches 
can be found in Section 2. 

The impact on the space of SUSY models has traditionally been presented in rather constrained frame¬ 
works, which have particular limitations when considering large numbers of analyses. One frequently 
used strategy for interpretation is in terms of models motivated by a particular mechanism of SUSY- 
breaking, for example via gravitational or gauge interactions. While such models can have theoretically 
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appealing features, they assert relationships between SUSY-breaking parameters that may not be realised 
in nature, and they sample only a small part of the parameter space of the MSSM. SUSY searches at the 
LHC have also been interpreted using ‘simplified’ models. Such models attempt to capture the behaviour 
of a small number of kinematically accessible sparticles, assuming all others play no role. The simplest 
case corresponds to one specific SUSY production process wifh a fixed decay chain. Such models provide 
insighf info fhe experimenfal consfrainfs on fhe individual sparficle and decay mode, buf fail fo capfure 
fhe complex effecls fhaf can resulf from large numbers of compefing production and decay processes. 

The MSSM has over a hundred paramefers fhaf describe fhe paffem of sparficle masses and fheir decays. 
This paramefer space is foo large fo be scanned exhausfively and compared fo ATLAS dafa. By apply¬ 
ing a series of assumpfions motivated by eifher experimenfal consfrainfs or general fealures of possible 
SUSY breaking mechanisms, fhe number of paramefers can be reduced fo 19. This is known as fhe phe¬ 
nomenological MSSM (pMSSM) [32-34]. This model is assumed fo conserve R parify,^ which ensures 
fhaf sparticles are produced in pairs and fhe lighfesf supersymmefric parficle (LSP) is sfable. The para¬ 
meters are assumed fo be real so fhaf new CP violation does nol occur in fhe sparficle secfor. Paramefers 
fhaf would give rise fo additional flavour-changing neufral currenfs are absenf. The LSP provides a dark 
mailer candidale if colourless and eleclrically neufral. In Ihis paper, fhe LSP is required fo be fhe lighfesf 
neulralino. Ifs producfion al fhe LHC gives rise fo missing Iransverse momenlum (whose magnilude is 
denoled which is required by mosl of fhe ATLAS searches considered in Ihis paper. There is no 

Iheorelical upper bound on fhe parameters characlerising fhe sparficle masses. However, since fhe ex- 
perimenfs have no sensifivily fo sparticles wifh very large masses, fhe following addilional restriclion is 
applied before a specific sel of paramefers is considered: all sparficle masses musl be less fhan 4 TeV. A 
specific sel of fhe 19 paramefers is referred fo as a model poinf in parameter space; 310,327 such model 
poinls, each consisfenl wifh a range of previous experimenfal resulls, are considered. More delails aboul 
fhe selection of pMSSM poinfs can be found in Section 3. 

Several groups have advocated fhe use of fhe pMSSM for inlerprelafion of LHC resulls [33-52]. Mosl 
of Ihese sludies use eslimafed experimenfal efficiencies and accepfances for pMSSM poinls, and com¬ 
pare Ihem fo fhe model-independenl limils from a selecfion of LHC searches fo conslrain fhe pMSSM 
parameter space. Previous ATLAS analyses have also used fhe pMSSM for inlerprelafion of individual 
searches [53-57] by fixing mosl of fhe paramefers, and varying jusl Iwo or Ihree: Ihey Iherefore explore 
only a small part of fhe parameter space. This paper makes full use of fhe ATLAS experimenfal simula- 
lion, reconslruclion and analysis fools. If represenls fhe mosl comprehensive assessmenl of fhe ATLAS 
consfrainfs on supersymmefry models fo dale. 

The paper is organised as follows. The relevanf ATLAS Run 1 analyses are summarised in Section 2. A 
descripfion of fhe pMSSM paramefer space can be found in Section 3, along wifh fhe direcf and indirecl 
consfrainfs applied prior fo fhe generalion of fhe 310,327 model poinls. Monle Carlo simulation of Ihose 
model poinls is described in Secfion 4. The effecl of fhe ATLAS searches on fhis pMSSM space is 
described in Section 5. Discussion and conclusions can be found in Secfion 6. 


2. ATLAS searches 

A folal of 22 dislincl ATLAS analyses are considered, spanning a wide range of differenl search sfrafegies 
and final slates, as lisfed in Table 1. Each analysis has several signal regions — for example fhe analysis 

* R = (_i)3®+i+2S where B is baryon number, L is lepton number and S is spin. 
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Analysis 

Ref. 

Category 

0-leplon -r 2-6 jels -t F™"" 

[58] 


0-leplon -t 7-10 jels - 1 - F™*® 

[59] 


1 -leplon + jels -i- 

[60] 


TijH) + jels -t 

[61] 

Inclusive 

SS/3-lepfons -t jels -t 

[62] 


0/1-leplon - 1 - 3?>-jels -i- 

[63] 


Monojel 

[64] 


0-leplon slop 

[65] 


1 -leplon slop 

[56] 


2-leplons slop 

[66] 


Monojel slop 

[67] 

Third generation 

Slop wilh Z boson 

[68] 


2/?-jels - 1 - 

[69] 


tb+E™^^, slop 

[57] 


ih 

[70] 


2-leplons 

[54] 


2-t 

[55] 

Eleclroweak 

3-leptons 

[53] 


4-leptons 

[71] 


Disappearing Track 

[72] 


Long-lived particle 

[73, 74] 

Olher 

HI A T+T- 

[75] 



Table 1: The 22 different ATLAS searches considered in this summary paper. The term ‘lepton’ {{) refers specific¬ 
ally to and states, except in the cases of the electroweak 3-leptons and 4-leptons analyses where t leptons are 
also included. 


requiring events with zero isolated electrons and muons and a minimum of 2-6 jets in association with 
large [58] has 15 different signal regions, each with different requirements on kinematic parameters 
and/or multiplicities of jets. For each of the 22 analyses, most of the signal regions from the original 
analysis are considered. However, in some cases, for practical reasons it was necessary to leave out some 
specialised signal regions or more complex combined fits. This leads to a slight underestimate of the full 
reach of the search. In total, almost 200 distinct signal regions are considered. 

The analyses are classified info fhe four broad cafegories shown in Table 1. ‘Inclusive’ searches are fhose 
primarily fargefing decays, including cascade decays, initialed by producfion of squarks of fhe firsl Iwo 
generalions or gluinos. ‘Third-generalion’ searches are fhose largeled particularly af fhe producfion of lop 
and bolfom squarks, known as fhe slop (?) and sbolfom (b). ‘Elecfroweak’ searches include fhose for direcl 
producfion of eleclroweakinos and slepfons. Since each search involves multiple signal regions, and since 
differenl SUSY producfion and decay processes can confribule lo each of fhose, Ihis calegorisafion can 
only be considered lo be a rough guide when interpreting fhe lype of sparticles lo which fhe analysis mighl 
show sensilivily. ‘Olher’ searches are fhose for heavy, long-lived parlicles (which are only considered for 
a small subsel of fhe model poinls) and fhe search for heavy Higgs bosons. The delails of fhe analyses can 
be found in fhe corresponding papers (listed in Table 1) and a brief summary for each is given below. 
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In what follows the term ‘lepton’ {£) is used to refer specifically to the charged leptons e- and jj.- of the 
first two generations. Where r- leptons are also included — for the 3-leptons and 4-leptons electroweak 
searches — this is indicated explicitly. 


2.1. Inclusive searches 

The inclusive searches are designed to be sensitive to prompt decays of squarks, particularly those of the 
first two generations, and gluinos. Strongly interacting sparticles may decay directly to the LSP, via the 
decay q ^ q +x^ for the squark and via g q + q for the gluino. Alternatively, cascade decays may 
also occur involving one or more additional sparticles yielding final states with additional jets, large 
and possibly leptons, including r leptons. The ATLAS searches targeting these final stales are classified 
according to the different dominant signal signatures, as follows. 

The 0-lepton -i- 2-6 jets -i- analysis [58] has wide-ranging sensitivity to strongly interacting sparticle 
production. It vetoes events with leptons in order to suppress the background from W boson and tt 
decays. Depending on the signal region, final slates with a minimum jet requirement of 2-6 jets with 
large transverse momenta (denoted px in the following) are considered, each in association with large 
Signal regions with small numbers of jets provide sensitivity to the direct production and decay 
of squarks, while those with higher jet multiplicities are sensitive to the production and direct decay of 
gluinos, and various cascade decays. The original analysis has two signal regions specifically largeling 
hadronic decays of high-px W bosons. These signal regions are not considered in this paper. 

The 0-lepton -i- 7-10 jets -i- analysis [59] selects events with significant and with jet multipli¬ 
cities ranging from 7 to 10 or more, depending on the signal region. It was designed to target, amongst 
others, models where each gluino of a produced pair decays through a (possibly virtual) top squark to 
t + t + x^. The four top quarks produced generally lead to large jet multiplicities in the final stale. This 
search also has sensitivity to other models in which cascade decays generate large numbers of jets. It 
has a looser requirement on the than the 2-6 jet analysis described in the previous paragraph, be¬ 
cause the many possible intermediate stages of the cascade decay tend to reduce the Unlike the 

original analysis, where the disjoint signal regions could be statistically combined to improve sensitivity, 
in this paper signal regions are considered individually when determining whether or not a model point is 
excluded. 

The 1-lepton -i- jets -i- L"™®* analysis [60] explicitly requires one isolated lepton, several jets and high 
in the selection. Two sets of signal regions are used from this analysis - one set with relatively high 
px leptons (sensitive to SUSY scenarios with larger mass splittings between the produced sparticle and 
the LSP) and another set using low-px leptons (sensitive to smaller mass splittings). The signal regions 
requiring two leptons are not considered. Overall, it is sensitive to decay chains where leptons can be 
produced through the cascade decay of squarks and gluinos. 

The t{tI£) + jets -i- search [61] targets final states arising from cascades producing hadronically 
decaying r leptons - with signal regions requiring either one or two r leptons, and including large 
jets and either exactly zero or one additional light lepton. This search can be sensitive to long decay 
chains in models with light staus. 

Cascade decays of squark and gluino pairs can also lead to final states with multiple leptons, or with two 
leptons of the same electric charge, known as same-sign (SS) leptons. Those final stales are addressed 
by the SS/3-leptons -i-jets -i- L"™®* analysis [62], which requires multiple jets in the final slate, and either 
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two SS leptons - with or without jets containing Zj-hadrons (Zj-jets) in the final state - or at least three 
leptons. 

For models where many fj-jets are expected, the specially designed 0/1-lepton -i- 3/j-jets -i- ana¬ 
lysis [63] is sensitive. This analysis is designed around the definition of two sets of signal regions: one 
set with no isolated leptons, and another with at least one isolated lepton. 

The Monojet analysis [64] selects events where the leading jet’s px is as large as 50% of the and 
there is large F"™®® and no leptons. The single jet can originate from initial-state QCD radiation (ISR), 
providing sensitivity to collisions in which no decay products from sparticle decays are observed. This 
can occur either for direct pair production of invisible LSPs, or if the produced sparticles are only a little 
heavier than the LSP (up to a few GeV) and their decays therefore produce SM particles of too low an 
energy to be detected in the other searches. 


2.2. Third-generation searches 

This set of analyses is focused on searches for direct production of third generation squarks. Their masses 
are generally expected to be at the TeV scale or below if the Higgs boson(s) are to be protected from large 
unnatural loop corrections. The decay of t and b squarks also leads to distinctive experimental signatures, 
typically involving the production of t- or /j-quarks in association with large 

The 0-lepton stop search [65] is optimised for the direct production of pairs of top squarks decaying 
directly to a top quark and neutralino, leading to an all-hadronic final sfafe wifh af leasf fwo /?-jefs and 
large S™®®. Mosf of fhe signal regions rely on variables relafed fo fhe reconsfrucfed fop quarks presenf in 
fhe final sfafe and on lepfon vefoes, buf fhere are also signal regions fhaf fargef fhe case where one of fhe 
pair of fop squarks decays fo a fop quark and a neufralino and fhe ofher decays fo a boffom quark and a 
chargino. 

The 0-lepfon sfop search is complemenfed by fhe 1-lepfon sfop search [56], in which all of fhe signal 
regions are characferised by exacfly one isolafed lepfon, af leasf fwo jefs and large The presence 

of a b-]e.t is used in bofh signal regions fargefing t\ —> bx\ and fhose fargefing t\ —> tx\, while fhe laffer 
also use variables relafed fo reconsfrucfed fop quarks. The dedicafed signal regions fargefing fop squark 
decays wifh soff lepfons or boosted top quarks are not included in this paper. 

The 2-leptons stop search [66] is designed for final sfafes confaining fwo isolafed lepfons and large 
primarily fargefing fop squarks decaying fhrough an infermediafe chargino. Only fhe so-called “lepfonic 
WIT 2 ” signal regions, fargefing charginos decaying through on-shell W bosons, are included in this pa¬ 
per. 

The Monojet stop analysis [67] looks for final sfafes characferised by large af leasf one high-px jet 
(vetoing events with more than three jets), and no leptons. The signal regions of this search were designed 
in the context of a search for top squarks, each decaying into an undetected charm quark and a neutralino 
which is relevant for cases where the LSP mass is close to the top squark mass. The requirement 
is less stringent than that of the Monojet analysis described in Section 2.1. The signal regions targeting 
reconstructed and tagged charm jets are not included in this paper. 

The search for top squarks with a Z boson in the final sfafe [68] is mofivafed by fhe decay of h —> tiZ, 
which can produce many lepfons in fhe final stale. The lepfons are required fo form a pair wifh a mass 
consislenf wifh fhe Z boson, wifh af leasf one b-]e.t and large 
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The lb-]eX& + analysis [69] searches for SUSY scenarios that produce events containing exactly two 
Z 7 -jets, significant and no isolated leptons, for example those coming from decays of bottom squarks 
to a fj-quark and the LSP and from top squark to a fj-quark and chargino. 

The tb+E™^^, stop analysis [57] was designed for a mixed scenario: direct production of pairs of top or 
bottom squarks each decaying (with various branching ratios) to neutralinos or charginos, and yielding 
final states consisting of a top quark, bottom quark and large 

2.3. Electroweak searches 

This section details the analyses considered in this paper which target sparticles produced via electroweak 
interactions. This includes the production of pairs of sleptons or electroweakinos which typically decay 
into final states containing several high-px leptons and significant E™^^. 

The lepton plus Higgs boson {ih) analysis [70] is designed to search for direct pair production of a 
chargino and a neutralino, which decay to final states with large S™*®, an isolated lepton, and a Higgs 
boson h which is identified by requiring either two fj-jets, or two photons, or a second lepton with the 
same electric charge (targeting h WW decays). Only the signal regions for Higgs boson decays to 
bottom quarks are considered in this paper. 

The 2-leptons analysis [54] targets electroweak production of charginos and/or neutralinos, or sleptons 
in events with exactly two leptons, large E^™^^ and, for some signal regions, two or more jets in the final 
state. 

A complementary search targeting the third-generation leptons is the 2 -t analysis [55] searching for 
SUSY in events with at least two hadronically decaying t leptons, large and a jet veto. 

The 3-leptons analysis [53] is a search for the direct production of charginos and neutralinos in final 
states with three leptons — which here may include up to two hadronically decaying t leptons — and 
large which can come through the decays via sneutrinos, sleptons or W, Z or Higgs bosons. 

The 4-leptons analysis [71] looks for SUSY in events with four or more leptons, of which at least two must 
be electrons or muons. The leptons may also include hadronically decaying r leptons in this case. Such 
high lepton multiplicity final states can occur if a degenerate T 3 P^ir is produced which subsequently 
decay via sleptons, staus or Z bosons to^® and many leptons. 

The Disappearing Track analysis [72] is motivated by scenarios with a wino-like LSP in which the charged 
wino is typically only ~ 160MeV heavier than the LSP. In such models the can have decay lengths of 
order a few tens of centimetres before it decays to a)f j and a charged pion. The low-momentum pion track 
is typically not reconstructed, so the distinctive signature is that of the high-px chargino track apparently 
disappearing within the detector volume. 

2.4. Other searches 

The long-lived particle searches [73, 74] are designed to detect heavy long-lived particles by measuring 
their speed using the time-of-flight to the calorimeters and muon detectors and jiy (where y is the 
relativistic Lorentz factor) from the specific ionisation energy loss in the pixel detector. Only the direct 
production of pairs of long-lived top or bottom squarks, gluinos, staus or charginos are considered in this 
paper. The search using 7 TeV data from 2011 [73] considered sparticles as light as 200 GeV, whereas in 


most cases the later analysis [74] only considered sparticles above 400 GeV. Both searches are therefore 
included for maximal sensitivity. 

The HjA —> T*T~ search [75] is designed to detect the heavy, neutral Higgs bosons predicted in the 
MSSM if they decay to r-pairs. 


3. pMSSM points and indirect constraints 

An overview of the method of selecting signal model points is provided in this section. A summary of the 
scan of the pMSSM parameter space, the software employed, and the constraints applied to determine the 
final selection are described. 


3.1. pMSSM points generation 

The model set is generated by selecting model points within the pMSSM using methods similar to those 
described in Ref. [34], but with several important changes. The modifications are made after taking 
into account new experimental results, updated calculational tools, knowledge gained from the study 
described in Ref. [34], and the improved capabilities of the ATLAS simulations. The full details of the 
method by which the model points are selected, including the sampling procedure, the codes employed, 
and the constraints applied to determine the final selection of ‘surviving’ model points may be found in 
the Appendices. 

The model points are selected after making the following assumptions about the MSSM. They are motiv¬ 
ated both by constraints from experimental observations and a desire for theoretical simplicity: 

(i) R-parity is exactly conserved. 

(ii) The soft parameters are real, so that no new sources of CP violation exist beyond that present in the 

CKM matrix. 

(ill) Minimal Flavour Violation [76] is imposed at the electroweak scale. 

(in) The first two generations of squarks and sleptons with the same quantum numbers are mass degener¬ 
ate, and their Yukawa couplings are too small to affect sparticle production or precision observables. 

(n) The LSP is the lightest neutralino. 

This approach remains agnostic about the presence of non-minimal particle content at higher scales, the 
mechanism of SUSY breaking, and the unification of sparticle masses. Assumptions ii-iv are motivated 
by the necessity of imposing some organising principle on SUSY flavour-violating parameters to allow 
TeV-scale masses for the squarks and sleptons. Combining assumptions i-v reduces the large MSSM 
parameter space to the 19-dimensional subspace considered here. The parameters and the ranges used 
to sample them are listed in Table 2. The 4 TeV upper bound on most of the mass parameters is chosen 
to make all states kinematically accessible at the LHC. As might be expected, decreasing the value of 
this upper limit restricts the space, resulting in an increase in the apparent fraction of the pMSSM space 
to which ATLAS analyses are sensitive. Further increasing any physical mass above 4 TeV has little 
effect on the LHC phenomenology in most cases. An exception is when decays proceed via virtual heavy 
sparticles, when increasing that sparticle mass would lead to further suppression of those decays. A larger 


9 


Parameter 

Min value 

Max value 

Note 


90GeV 

4TeV 

Left-handed slepton (first two gens.) mass 

rrieii^ nie^) 

90GeV 

4TeV 

Right-handed slepton (first two gens.) mass 

m-L, 

90GeV 

4TeV 

Left-handed stau doublet mass 

m-ci 

90GeV 

4TeV 

Right-handed stau mass 

mQ^{= niQ^) 

200 GeV 

4TeV 

Left-handed squark (first two gens.) mass 

mui(= mu2) 

200 GeV 

4TeV 

Right-handed up-type squark (first two gens.) mass 

>ndi(= 

200 GeV 

4TeV 

Right-handed down-type squark (first two gens.) mass 


100 GeV 

4TeV 

Left-handed squark (third gen.) mass 

mui 

100 GeV 

4TeV 

Right-handed top squark mass 


100 GeV 

4TeV 

Right-handed bottom squark mass 

IMil 

OGeV 

4TeV 

Bino mass parameter 

IM2I 

70 GeV 

4TeV 

Wino mass parameter 

lyu| 

80 GeV 

4TeV 

Bilinear Higgs mass parameter 

M3 

200 GeV 

4TeV 

Gluino mass parameter 

|Ad 

OGeV 

8 TeV 

Trilinear top coupling 

|Afo| 

OGeV 

4TeV 

Trilinear bottom coupling 

lArl 

OGeV 

4TeV 

Trilinear r lepton coupling 

Ma 

100 GeV 

4TeV 

Pseudoscalar Higgs boson mass 

tanyS 

1 

60 

Ratio of the Higgs vacuum expectation values 


Table 2: Scan ranges used for each of the 19 pMSSM parameters. Where the parameter is written with a modulus 
sign both the positive and negative values are permitted. In the above, “gen(s)” refers to generation(s). 


range is permitted for |A,|, a parameter which affects loop corrections to the mass of the the Higgs boson. 
The larger range increases the fraction of model points having the mass of the lightest Higgs boson close 
to the measured value. 

Given the large dimensionality of the pMSSM, a grid sampling technique at regular intervals is imprac¬ 
tical. The space is therefore sampled by choosing random values for each parameter. It should be noted 
that in many cases only some of the parameters are relevant for a given observable, in which case the scan 
is effectively more comprehensive within the subspace of relevant parameters. The value of each para¬ 
meter is chosen from a flat probability distribution, with lower and upper bounds given in Table 2. The 
lower and upper limits of the parameter ranges are chosen to avoid experimental constraints and to give a 
high density of model points with masses at scales accessible by the LHC experiments, respectively. 

Condition iv imposes the constraints that the soft mass terms for the second generation are equal to those 
in the first, as shown in Table 2. This means, for example, that ml and cl have the same soft mass term 
in the Lagrangian so that their physical masses are very close. Furthermore the scalar partners of the 
left-handed fermions, such as and have the same soft mass due to S U{2)i^ invariance, but D-terms 
related to electroweak symmetry breaking split their mass-squared values by 0{m^). 

Once each of the 19 parameters has been chosen, a variety of publicly available software packages are 
used to calculate the properties of each model point, as described in Appendix A. In some cases the 
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Parameter 

Minimum value 

Maximum value 

Ap 

-0.0005 

0.0017 

> 

1 

-17.7 X 10“^° 

43.8 X 10-1° 

BR(t7 — > 5y) 

2.69 X 10-^ 

3.87 X 10-4 

BR(R, ^ p+p-) 

1.6 X 10-'^ 

4.2 X 10-° 

BR(B+ ^ T+v^) 

66 X lO”'^ 

161 X 10“° 

A 1 

— 

0.1208 

r invisible(SUS Y) (^) 

— 

2MeV 

Masses of charged sparticles 

100 GeV 

— 

m(xf) 

103 GeV 

— 

m{Ui2, ^1,2. Cl,2, 5i,2) 

200 GeV 

— 

m(h) 

124 GeV 

128 GeV 


Table 3: Constraints on acceptable pMSSM points from considerations of precision electroweak and flavour results, 
dark matter relic density, and other collider measurements. A long dash (—) indicates that no requirement is made. 
Further details may be found in the text. 


software is modified to produce accurate results for the wide range of models found in the pMSSM scan. 
The sparticle decays are calculated, again using a variety of codes and analytical techniques, as described 
in Appendix B. 


3.2. pMSSM point selection 

Acceptable model points are furthermore required to have consistent electroweak symmetry breaking, 
a scalar potential that does not break colour or electric charge, and all particles’ mass-squared values 
must be positive. Model points with theoretical pathologies, described in more detail in Appendix C, 
are discarded. Further experimental constraints, shown in Table 3, which indirectly affect the parameter 
space are applied and described below. 


3.2.1. Precision electroweak and flavour constraints 

Unless specified otherwise, the relevant observables are calculated using micrOMEGAs 3.5.5 [77, 78]. 
The constraint on the electroweak parameter Ap uses the limit on AT (the parameter describing the radi¬ 
ative corrections to the total Z boson coupling strength, the effective weak mixing angle, and the W boson 
mass) in Ref. [79] and Ap = a AT with a = 1/128. The allowed branching ratio (BR) of b ^ sy is the 
union of the two standard deviation (2cr) intervals around the theoretical prediction and the experimental 
measurement from Ref. [80]. For the branching ratio of Bs the value calculated by micrOMEGAs 

is scaled by 1/(1 - 0.088) as proposed in Ref. [81] for comparison with experiment. The scaled value 
is required to lie within the 2cr interval around the combined result from the LHCb and CMS Collab¬ 
orations [82]. The 2cr theoretical prediction for the SM (3.20 to 4.12) x 10“^ lies within this interval. 
The branching ratio is calculated using Ref. [83], which includes tanyS-enhanced corrections. 
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The allowed range is the union of the 2cr intervals around the experimental results [84-87] and the SM 
prediction [88]. Finally, for the SUSY contribution to the anomalous magnetic moment of the muon, 
A(g - 2)^, a very large range is allowed. This range in A{g - 2)^ is the union of the 3cr intervals around 
the SM value, (0.0 ± 5.9) x 10“^** from combining [89] the results of Refs. [90-92] and the experimental 
measurement [93, 94] corrected to an updated value of the muon to proton magnetic moment ratio [95, 
96] giving A{g - 2)p exp = (24.9 + 6.3) x 10~^°. Three-sigma intervals are used to obtain a continuous 
range from the union. 


3.2.2. Dark matter constraints 

Since R-parity is assumed to be exactly conserved, the LSP is stable and as a consequence has a non¬ 
zero cosmological abundance. It is assumed that the LSP abundance is determined thermally and is 
not diluted by other processes e.g. late-time entropy addition. No assumption is made about whether 
the LSP is the sole constituent of dark matter. As a result, the total cold dark matter energy density is 
used as an upper limit on the LSP abundance. The limit is based on the latest combined measurement 
from the Planck Collaboration of Q.cDMh^ = 0.1188 ± 0.0010 (Table 4 of Ref. [97]).^ The upper limit 
is set to the observed central value plus double the experimental uncertaintity. The limit on the spin- 
independent cross-section is that for the interaction of a neutralino with a nucleus derived by the LUX 
experiment [98]. In the case of the LSP mass versus proton spin-dependent cross-section plane the limit is 
from the COUPP Collaboration [99], while in the LSP mass versus neutron spin-dependent cross-section 
plane, the XENONIOO Collaboration [100] limit is applied. MicrOMEGAs 3.5.5 is used to calculate 
the neutralino-nucleon cross-sections. These are scaled down by the ratio of the expected relic density 
from the LSP to the observed relic density to obtain the effective dark matter cross-sections, assuming 
the remaining non-LSP dark matter is invisible to the direct detection experiments. When accepting or 
rejecting models, the calculated value is allowed to be up to a factor of four higher than the limits obtained 
by the experiments, to account for nucleon form-factor uncertainties [33]. 


3.2.3. Collider constraints 

Finally, constraints from LEP and from the measurement of the Higgs boson mass at the EHC are applied. 
To ensure consistency with EEP, model points are discarded if their additional contribution to the invisible 
width of the Z boson is above 2 MeV [101], or where any charged sparticle is lighter than 100GeV. Eor 
charginos, the bound is increased to 103 GeV, provided that all sneutrinos are heavier than 160 GeV 
and the mass splitting between the chargino and the ESP is at least 2 GeV. This constraint comes from 
the combined EEP search [102]. Eirst- and second-generation squarks are required to be heavier than 
200 GeV, although this has only a very small effect given the scan range and the assumption of negligible 
first- and second-generation squark mixing. The lightest Higgs boson mass, as calculated by FeynHiggs 
2.10.0 [103, 104], is required to be in the range 124 to 128GeV. This range is set around the central 
value of the Higgs mass at the time of generation, 126 GeV, and with a 2 GeV window that mainly reflects 
the typical theoretical uncertainty of the FeynHiggs calculation. The results are found not to depend on 
the exact value of the Higgs mass within this interval (as shown later in Section 5.4). 


^ It should be noted that in the context of the dark matter relic density the symbol h corresponds to the normalised Hubble 
constant, rather than the Higgs boson. 
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LSP type 

Definition 

Sampled 

Simulated 

Number Fraction 

Weight 

‘Bino-like’ 


480 X 10^ 

103,410 

35% 

1/24 

‘Wino-like’ 


)20 X10 ^{ 

80,233 

26% 

1 

‘Higgsino-like’ 


126,684 

39% 

1 

Total 


500 X 10^ 

310,327 




Table 4; Categorisation of the 310,327 model points by the type of the LSP (assumed to be the^j) according to 
the neutralino mixing matrix parameters Nij, where the hrst index indicates the neutralino mass eigenstate and the 
second indicates its nature in the lexicographical order {B, W,H\,H 2 )- For example, N\;i is the amplitude for the 
LSP to be W. The hnal two columns indicate the fraction of model points in that category that are sampled, and 
their weighted fraction after importance sampling. 


3.2.4. Importance sampling by LSP type 

Since low-mass SUSY models typically over-produce dark matter, the relic density constraint in Table 3 
sculpts the distribution of the allowed model points. The constraint depends strongly on the nature of the 
LSP. Except where particularly effective neutralino annihilation mechanisms are available, model points 
with a bino-like LSP generally tend to produce too much dark matter [105], meaning that such models 
are infrequently sampled and accepted in a random scan employing flat priors. The model points are 
therefore partitioned into three categories, bino-like, wino-like and Higgsino-like. The categorisation is 
made according to the dominant contribution to the LSP within the neutralino mixing matrix Nij as shown 
in Table 4. Model points are therefore selected, by importance sampling, in such a way that approximately 
equal numbers are obtained for each LSP type. In total 500 million model points are sampled randomly 
within the ranges listed in Table 2. From the first 20 million sampled, 206,917 model points had a wino- 
like or Higgsino-like LSP and satisfied all of the constraints of Table 3. To obtain a sufficiently high 
number with bino-like LSP, the remaining 480 million model points are used to find the 103,410 which 
had a bino-like LSP and satisfied the Table 3 constraints. Generally models have a LSP dominated by one 
particular type, with over 87% of models having a LSP which is at least 90% pure. The phenomenology 
of each LSP type can be explored separately due to the large number of model points in each category. 
In the following plots, where all LSP types are shown together the contribution from each LSP type is 
scaled according to the weights shown in Table 4. 


3.3. Properties of model points (before applying ATLAS constraints) 

The distributions of the gluino and LSP masses for the model points satisfying the constraints from Table 3 
are shown in Figure 1, separately for models with a bino-, wino- or Higgsino-like LSP. Light gluinos are 
more common among model points with bino-like LSP. Dark matter for model points with bino-like LSP 
is typically over-produced, so the presence of a gluino state close in mass to the LSP enables them to 
act as coannihilators with the dark matter in the early universe, reducing the relic density to a level that 
satisfies the constraint on The neutralino mass distribution for the bino-like LSP model points 

shows a sharp concentration of model points with m(x^) < 100 GeV. This concentration corresponds to 
model points in which the dark matter relic density constraint can be satisfied due to enhanced neutralino 
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Figure 1: Distributions of the gluino and LSP masses, shown separately for models with a bino-like (dotted red), 
wino-like (dashed blue) or Higgsino-like (solid green) LSP. The constraints listed in Table 3 have been applied, but 
not the constraints from the ATLAS searches. The distributions are normalised to unit area. The inset in the plot on 
the right shows in more detail the region of low neutralino mass for the models with bino-like LSP. 



Figure 2: Normalised distribution of sparticle masses for the lightest top (solid red), bottom (dashed blue), first- or 
second-generation squark (dashed green) and stau (dot-dashed black) for all LSP types combined. The constraints 
listed in Table 3 have been applied, but not the constraints from the ATLAS searches. 
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pair annihilation via the Z or Higgs boson poles. A plot of the low-mass region for bino-LSP model 
points (shown in the inset to Figure 1(a)) confirms this interpretation, showing two individual peaks 
corresponding to the separate Z and Higgs poles. 

The identity of the next-to-lightest supersymmetric particle (NLSP) can be important for the phenomeno¬ 
logy of a model point, and is strongly influenced by the LSP type. The NLSP is nearly always a chargino 
or neutralino for wino-like or Higgsino-like LSPs, as is expected given the small splittings between the 
different components of the wino or Higgsino multiplets. In particular, over 99% of wino-like LSP model 
points have a chargino NLSP. The small mass difference between the and the can lead to long xf 
lifetimes for wino-like and Higgsino-like LSPs, which result in the types of disappearing tracks searched 
for in the analysis described in Ref. [72]. 

Bino-like LSP model points, by contrast, exhibit a much wider range of NLSP types, the distribution 
of which is determined by their effect on the LSP annihilation rate. Even for models with a bino-like 
LSP, charginos and neutralinos are still the most common NLSP type, since LSP-Higgsino mixing is 
important for many of the possible annihilation mechanisms. The remaining NLSP types are generally 
ordered by their effectiveness as coannihilators, with coloured NLSPs being the most prevalent and neutral 
NLSPs (sneutrinos) being the least so. The NLSP-LSP mass splitting for bino-like LSP model points is 
frequently less than 50 GeV, supporting the assertion that light sparticles beyond the LSP are typically 
required to avoid over-production of bino-like LSPs. 

In cases where a squark, gluino or slepton is almost mass degenerate with the LSP, the small available 
phase-space in the decay can lead to sparticles that are long-lived or even stable on the time scale for 
traversal of the ATLAS detector. In total 3,427 model points contain squarks, gluinos or sleptons with 
CT > 1 mm. The sensitivity of the SUSY searches targeting prompt decays (the first three categories in 
Table 1 , other than the Disappearing Track analysis) to these model points is reduced and therefore only 
the long-lived particle searches were considered for these model points. 

The distributions of the mass of several other sparticles of phenomenological interest at the LHC can be 
found in Figure 2. Even before ATEAS analyses are considered, the model set is depleted in light top 
squarks, and to a slightly lesser extent in light bottom squarks, by the requirement (in Table 3) that the 
lightest Higgs boson mass be close to the experimentally observed value. 


4. Signal simulation and evaluation of searches 

For each of the 310,327 model points passing the preselection described in Table 3 it has to be determined 
which, if any, of the ATEAS searches are sensitive to it and whether it can be excluded or not. Simulating 
and running the full set of ATEAS analyses on these would be extremely time and resource consuming. 
Each model point is therefore evaluated in three steps as described below. A special procedure is applied 
to evaluate the searches for long-lived particles and heavy Higgs bosons as described in the following 
sub-sections. 


4.1. Supersymmetry signals 

First, each model point is categorised based on its production cross-sections for SUSY particles, as cal¬ 
culated using Prospino 2.1 [106-110]. The production processes are split into three separate groups: 
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Production mode 

Minimum cross 

section [fb] 

Fraction of models generated 

Bino LSP 

Wino LSP 

Higgsino LSP 

Strong 

0.25 

82.5% 

74.9% 

76.7% 

Mixed 

0.25 

52.6% 

42.1% 

13.9% 

Electroweak 

7.5 

38.3% 

72.5% 

75.0% 

Slepton pair 

0.75 

9.6% 

7.9% 

9.5% 


Table 5; Minimum cross-sections required to do particle-level event generation for the four different production 
modes and the fraction of the models above this cross-section for each LSP type. 


strong production, electroweak production (encompassing electroweakino and slepton pair production) 
and finally mixed production (e.g. of an electroweakino in association with a squark or gluino). Model 
points with cross sections for any of those processes larger than the minimum values in Table 5 are sub¬ 
sequently retained, and any such processes are investigated in more detail as described below. For strong 
and mixed production, the minimum cross-section corresponds to just five signal events produced in the 
full ^/s = S TeV dataset. Sensitivity to such small cross-sections occurs only for model points with a 
very high fraction of events with four leptons in the final state, for which the 4-leptons analysis has a high 
acceptance (up to 50%). Production of electroweakinos is most effectively observed by using decays to 
leptons, which are often suppressed by the leptonic branching fractions of W, Z and h bosons, explaining 
the higher cross-sections limits. ATLAS searches have a greater sensitivity to low cross-section slepton 
pair production than to electroweakino production. Therefore, if the model point does not satisfy the 
higher cross-section criterion of the electroweak production group, a fourth group that allows for model 
points with lower cross-section slepton pair production is considered. Fewer than 10% of the model points 
have no process passing this selection and are thus not considered to be excludable. 

For each of the model points satisfying one or more of the production mode cross-section criteria, a large 
sample of events is generated using MadGraphS 1.5.12 [111] with the CTEQ 6L1 parton distribution 
functions (PDF) [1 12] and Pythia 6.427 [1 13] with the AUET2B set of parameters [1 14]. MadGraphS is 
used to generate the initial pair of sparticles and up to one additional parton in the matrix element, while 
Pythia is used for all sparticle decays and parton showering^ MLM matching [115] is used with up to 
one additional jet in the MadGraph matrix element, a MadGraph kT measure of 100 GeV, and a Pythia 
jet measure cut off of 120GeV. Both Tauola 1.20 [116] and Photos 2.15 [117] are enabled to handle 
the decays of r leptons and final-state radiation of photons, respectively. 

To reduce the amount of computationally expensive detector simulation and reconstruction that is re¬ 
quired, a Monte Carlo particle-level selection corresponding to each of the SUSY searches in the first 
three categories in Table 1 is used to process the generated events. In this step, inefficiencies from the 
detector-level reconstruction are parameterised using a single efficiency factor for each signal region, de¬ 
termined from previously simulated signal samples. Exceptions to this are made for the r reconstruction 
efficiency in the t(t/^) -i- jets -i- E™*® and 2-t searches, for which pj-dependent efficiencies are applied 
for each signal t. Similarly, the Disappearing Track search applies the reconstruction efficiency for de¬ 
caying charginos as a function of the distance from the centre of the ATLAS detector and the angular 
coordinates as published in Ref. [72]. The expected event yield in each signal region is calculated for 
each model point. For most analyses the categorisation is performed by directly comparing the expected 


^ Polarisation from the decay of the initial sparticles is not taken into account in this analysis. 
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signal yield A^sig to the model-independent 95% confidence level (CL) upper bound on the number of 
beyond-the-SM events in each signal region of that analysis. Model points are then partitioned into 
three categories, on the basis of that particle-level simulation, using criteria determined to be appropriate 
for each individual analysis. The first category comprised those already excluded at this stage on the 
basis that A^sig is sufficiently larger than A^^ax for at least one signal region of one analysis. The expected 
sensitivity of all other analyses to such model points is calculated using particle-level yields, and using 
average reconstruction efficiencies. The second category corresponds to those found not to be exclud¬ 
able, consisting of points with A^sig materially smaller than A^^ax for all analyses. The exact relationship 
between A^sig and A^^ax for tho categorisation is determined separately for each signal region and depends 
on the accuracy with which the particle-level evaluation reproduces the results of a full simulation. In 
total 35.9% (44.7%) of the model points fall in the first (second) category and are deemed to have been 
excluded (not excluded) at the 95% CL. The validity of this classification was confirmed using the full 
simulation and reconstruction procedure described below for approximately 5% of the model points in the 
first category for each analysis. A final category of model points - those with Nsig close to A^^ax (typically 
within a factor of a few) for the most sensitive analyses - are subject to more detailed investigation, as 
follows. 

For the 44,559 model points for which the overall exclusion is uncertain based on the particle-level sim¬ 
ulation described above, the final step is a fast, GEANT4-based [118] simulation using a parameterisation 
of the performance of the ATLAS electromagnetic and hadronic calorimeters [119] and full event recon¬ 
struction. The simulation includes a realistic description of multiple pp interactions per bunch crossing, 
and is corrected for identification efficiencies and resolution effects. For each such model point, signal 
events are generated corresponding to four times the integrated luminosity recorded (i.e. 81.2 fb“^). The 
simulation is limited to those production modes which could contribute to the analyses of interest. For 
these processes, the nominal cross-section and the uncertainty are taken from an envelope of cross-section 
predictions using different parton distributions and factorisation and renormalisation scales, as described 
in Ref. [120]. The addition of the resummation of soft gluon emission at next-to-leading-logarithm accur¬ 
acy (NLL) [106, 121-124] is performed in the case of strong production of sparticle pairs. 

The status of each previously inconclusive model point is then determined for each of the analyses using 
the same procedure [125] as used in the original analyses. In each analysis the signal region with the 
best expected sensitivity is identified and the “CLs method” [126] is used to determine if the model point 
is excluded or not at 95% CL. It should be noted that for the exclusion fits, the nominal signal cross- 
sections are used, without any theoretical uncertainties on the signal, except for the Monojet stop and 
Monojet analyses. These two analyses are particularly sensitive to the modelling of ISR as they rely on 
a high-pT ISR jet in their event selection. Therefore an additional 25% ISR signal uncertainty is applied 
in those cases, based on the observed variance in acceptance in signal samples with modified parameters 
for the ISR modelling [67]. For the 0-lepton -i- 7-10 jets -i- 0/1-lepton -i- 3A-jets -i- and 2- 

leptons stop analyses, it is not possible to apply the full combined fit procedure of the original analyses. 
Instead only the individual signal regions are considered, resulting in somewhat conservative modelling 
of the sensitivity for those three analyses. For the overall exclusion, no attempt is made to combine the 
individual analyses. Instead the analysis with the best expected exclusion is used for each model point to 
determine its status. 
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4.2. Long-lived particle search 

Model points with heavy long-lived particles require special treatment since such particles can traverse 
part or all of the ATLAS detector leaving rather distinct signatures. The dominant types of long-lived 
particles in the model points are the Xi and the X 2 when they are almost mass-degenerate with the LSP. 
The decay of such long-lived particles is included in the simulation procedure described above and model 
points with such long-lived particles are considered using the same procedure. 

Aside from the electroweakinos, 3,427 of the model points contain squarks, gluinos or sleptons with 
CT > 1 mm. These model points have not been simulated. Instead only the results from the Long-lived 
particle searches are used to constrain these model points. 

The long-lived particle searches in Refs. [73] and [74] provide limits on the production cross-section 
at 7 and 8 TeV, respectively, of bottom squarks, top squarks, gluinos, staus and charginos in the case 
where these live long enough to traverse the complete detector. Model points with bottom squarks, top 
squarks, gluinos, staus or charginos with a lifetime above 85 ns and production cross-sections exceeding 
the corresponding cross section limit from either the 7 or 8 TeV result are considered to be excluded. In 
all other cases - where the lifetime is shorter or the production cross-sections lower - the model point is 
considered not to be excluded. 

4.3. Heavy Higgs boson search 

Cross sections and branching ratios for heavy Higgs bosons are calculated for gluon fusion, or for pro¬ 
duction in association with Z7-quark(s) [127-129]. The high-ni^ category {niA > 200 GeV) of the ATLAS 
search [75] is used since this regime is relevant to all the model points in this study. It is assumed that 
Z7-quark associated production dominates and this calculation is performed for each model point in the 
pMSSM parameter space using the software SusHi 1.3.0 [103, 130-141]. The large value of ma in all 
the model points has the effect that the A and H bosons are nearly mass degenerate, so both must be 
simulated. The quantity o-{bbH) x BR(// —> rr) - 1 - cr{bbA) x BR(A —> tt) is calculated and is compared 
to the ATLAS 95% CL upper limits [75] for a scalar particle produced in association with f7-quark(s) and 
decaying to rr. For the overall exclusion, this heavy Higgs boson search is considered only if none of the 
SUSY searches are expected to exclude the model point at the 95% CL. 


5. ATLAS constraints from LHC Run 1 

5.1. Impact of ATLAS searches on sparticle masses 

The effect of the ATLAS search constraints are most easily presented as projections onto one-dimensional 
or two-dimensional subspaces of the full 19 parameter pMSSM space.^ The most relevant parameters 
onto which to project are typically the sparticle masses. Production cross-sections for sparticles decrease 
rapidly when their masses are increased. When those initial sparticles decay, the masses of other sparticles 
affect the types of visible decay products and their kinematics. The mass of the LSP is particularly 


^ A full list of model parameters, observables and which analyses, if any, are excluding each model, are available from the 
ATLAS Collaboration website [142]. 
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important since a decay to a high-mass LSP results in less energy being available for the observable 
decay products. 

The fraction of surviving model points in the projected space necessarily depends both on the prior dis¬ 
tribution of model points in the parameters that have been projected out, and on experimental constraints 
on sparticle masses other than those plotted. Thus, some care is needed in their interpretation. In par¬ 
ticular the fractions of model points excluded can depend, in some cases sensitively, on the non-collider 
constraints shown in Table 3, the choice of scan ranges shown in Table 2, and on the choice of a flat prior. 
Nevertheless, some general features of the impact of the ATLAS Run 1 searches are clear. 

The simplified-model limits shown on the plots throughout this section are the observed limits from the 
indicated analysis. In many cases there are several analyses interpreting their results in the same simplified 
models, and in this paper the observed limits from the most constraining analysis are shown. It should be 
noted that there is no minimum number of model points required in each bin. 

The results are shown first for squark masses (of the first two generations) and the gluino mass in Sec¬ 
tion 5.1.1, then for third-generation squark masses in Section 5.1.2, and for the electroweak sparticles in 
Section 5.1.3. A small subset of the model points contain long-lived squarks, gluinos or sleptons. These 
3,427 model points are treated separately in Section 5.1.4 as only the long-lived particle search is con¬ 
sidered for these model points. Section 5.1.5 describes the effect of a search for the decay of heavy neutral 
Higgs bosons to two t leptons. The complementarity between the different ATLAS searches is described 
in Section 5.1.6. 


5.1.1. Squarks and gluinos 

Figure 3(a) shows the fraction of model points excluded by the ATLAS searches as projected onto the 
two-dimensional space of the masses of the LSP and the gluino. As one would expect, light gluinos are 
robustly constrained by the ATLAS searches, whereas at larger gluino masses the fraction of model points 
excluded is reduced. 

It is instructive to compare these observed pMSSM exclusions to the observed limits previously presented 
for simplified low-scale models. Superimposed on Figure 3(a) is a line showing the 95% CL exclusion 
previously derived by the 0-lepton + 2-6 jets -i- analysis from a simplified model in which only 
the gluino and the LSP are kinematically accessible [58]. It can be seen that there is generally good 
congruence between the region excluded in the two different scenarios, demonstrating that the simplified 
model is successfully capturing the main pMSSM phenomenology in this case. Nevertheless, the pMSSM 
sensitivity does differ in detail from that of the simplified model. Not only is this because multiple 
analyses are considered here but also because of residual dependence on the masses of other sparticles. 
For example, the pMSSM permits the existence of additional particles with masses lying in between 
those of the gluino and the LSP, which can lead to cascade decays. In general such cascade decays can be 
expected to yield different jet pj and spectra, and perhaps additional leptons, which have different 
experimental acceptances. 

Close to the diagonal line where m{g) is only a little larger than m(x^) (Figure 3(a)) the simplified-model 
exclusion from Ref. [58] underestimates the ATLAS sensitivity. The g —> qqx^i decays in this near¬ 
degenerate region produce low-energy quarks which typically fail to meet the kinematic requirements 
on jets. That ATLAS does indeed show good sensitivity in this region is instead due to the monojet 
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Figure 3: Fraction of pMSSM points excluded by the combination of 8 TeV ATLAS searches in the (a) and 
(b) the mass planes. The colour scale indicates the fraction of pMSSM points excluded in each mass bin, with 
black squares indicating 100% of model points being excluded. The white regions indicate places where no model 
points were sampled which satisfied the constraints of Table 3. In both cases, the solid white lines overlaid are 
observed simplified-model limits from the 0-lepton + 2-6 jets + search [58] at 95% CL. In the case, the 
simplified-model limit is set assuming direct production of gluino pairs and that the squarks are decoupled, with 
gluino decaying to quarks and a neutralino, g —» q+q+x^y In the q-x\ plane, both lines are drawn assuming directly 
produced first/second-generation squark pairs, with each squark decaying to a quark and a neutralino, q ^ q +^[. 
The solid line corresponds to the case where all eight squarks from the first two generations are assumed to be 
degenerate. The dashed line has the squark production cross-section scaled down by a factor of four to emulate the 
effect of only two of those eight squarks being kinematically accessible. 


analyses [64, 67]. These analyses were designed to capture the recoil of LSPs (or other, slightly heavier, 
SUSY particles) against initial-state QCD radiation. 

Figure 3(b) shows a different projection, in this case to the mass of the LSP versus the mass of the lightest 
squark of the first two generations, ^l,r for q € {u, d, s, c), labelled here and in what follows as q. It can 
be observed that there is good sensitivity at low squark mass and no models with a squark mass below 
250 GeV are allowed by the ATLAS analyses. The solid line superimposed on Figure 3(b) shows the 
95% CL exclusion obtained previously [58] for a simplified model in which the only kinematically ac¬ 
cessible sparticles are the LSP and the eight squark states of the first two generations, where these squarks 
are all assumed to have the same mass. It can be seen that the region within the solid simplified-model 
exclusion curve is only partially excluded within the pMSSM. This is primarily because the pMSSM-19 
parameter space does not demand that the squarks be eight-fold degenerate, reducing the cross-section. 
There is a closer correspondence between the pMSSM sensitivity and that of an alternative simplified 
model (dashed line), in which the cross-section for direct (anti-)squark production has been reduced by a 
factor of four, to model the effect of only two of those eight squarks being mass degenerate.^ 

A noticeable excursion from the simplified-model lines, visible on both plots in Figure 3 is a horizontal 
band of sensitivity to pMSSM points for LSP masses less than about 200 GeV stretching up to large gluino 

^ Reference [58] emulates the effect of a single kinematically accessible squark by dividing the cross-section by a factor of 
eight rather than four. 
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Figure 4; Fraction of pMSSM points excluded in the g-q plane, where q represents the lightest squark from the 
first two generations. The overlaid line shows a limit for a simplified model from the 0-lepton + 2-6 jets + F™®* 
search [58] which assumes strong production of gluinos and eight-fold degenerate first- and second-generation 
squarks, with direct decays to quarks and massless neutralinos. The colour scale is as described in Figure 3. 


(or q) masses. Since such high-mass strongly interacting sparticles have small production cross-sections, 
one would not expect sensitivity to their production. Indeed these constraints are not the result of gluino 
or squark searches, but instead of searches for disappearing tracks from long-lived charginos. These 
long-lived chargino states are common for models with wino-like LSPs with mass splittings between 
the charged NLSP and the neutral LSP of less than about 200 MeV. The NLSP, when it decays inside 
the detector volume, produces an invisible LSP and a low-energy charged pion which itself often goes 
undetected. The search for such ‘disappearing’ charged-particle tracks is sensitive even in the absence 
of direct squark or gluino production, and hence sensitivity is observed for any mass of the strongly 
interacting sparticle. 

Figure 4(a) shows the sensitivity as projected onto the plane of the gluino and squark masses, where now 
the LSP mass may take any value. One observes near-total exclusion by ATLAS analyses of gluinos with 
masses less than about 700 GeV, with a high fraction of exclusion up to about 1.2 TeV, for all values of 
the lightest squark mass. Light squarks are also strongly constrained, although those constraints weaken 
as the gluino mass increases, due to suppression both of direct squark-pair production via t-channel gluino 
exchange and of associated production of ^ -i- 

The simplified model superimposed onto the squark-gluino plane is one that assumes an eight-fold de¬ 
generacy of squark masses in the first two generations and a massless LSP [58]. As one would expect, 
this simplified-model line lies close fo fhe upper edge of fhe pMSSM sensifivify, since fhe pMSSM per- 
mits non-degenerale squarks, and a non-zero LSP mass, bofh of which reduce sensifivify, by reducing fhe 
signal cross-secfion and experimenfal accepfance respectively. The reducfion in sensifivify caused by a 
non-zero LSP mass is more pronounced in fhe case of model poinfs wifh a bino-like LSP, Figure 4(b). 
These model poinfs offen have a small mass difference befween fhe squark and fhe LSP in order fo satisfy 
fhe dark matter relic consfrainf, as discussed earlier in Secfion 3.3. 
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Figure 5; Fraction of model points excluded in the planes of the masses of the left-handed and right-handed squarks 
(of the first two generations) versus the neutralino mass. Both simplified-model limit contours are taken from the 
0-lepton + 2-6 jets -H analysis (Figure 10(c) of Ref. [58]); however, for the left (right) handed squarks the 
assumption of four (two) degenerate squarks is emulated by dividing the cross-section for production of the eight 
degenerate squark states by a factor of two (four). The colour scale is as described in Figure 3. 
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The pMSSM also allows one to explore how the sensitivity to direct squark production depends on the 
nature of the lightest squark. Figure 5 shows that the search reach depends on whether the squarks are left- 
or right-handed, and whether they are up-type or down-type squarks. Thus the assumption of equivalent 
sensitivity for all such squarks, common in presentation of LHC SUSY searches, is a simplification that 
is not justified in fhe more general confexf of fhe pMSSM. The producfion of down squarks is suppressed 
relafive fo up squarks since fhere are fewer valence down quarks fhan up in fhe profon. The resulfs for 
lefl-handed squarks also differ from fhose for righf-handed squarks. The ml and squarks form a SU(2) 
doublef and so are degenerafe in mass up fo elecfroweak symmefry breaking effecfs. This means fhaf if 
fhe lighfesf squark is a ^l, there would be another squark with similar mass - a statement not usually true 
for ^R, which have no similar constraint. For a particular value of the lightest squark mass, the presence 
of a pair of left-handed squarks effectively increases the overall squark production cross-section, leading 
to larger apparent sensitivity. The improved sensitivity to left-handed squarks is further enhanced by the 
dominant squark decay modes. Right-handed squarks, which lack weak couplings and have small Yukawa 
couplings, have suppressed decays to the wino-like or Higgsino-like LSPs which dominate the model 
sample. Instead right-handed squarks generally cascade decay via other electroweakino states resulting 
in events with smaller a greater number of lower-px jets, and generally a smaller experimental 

sensitivity. 


5.1.2. Third-generation squarks 

The third-generation squarks are of particular phenomenological interest, since light top squarks (and to a 
lesser extent bottom squarks) are usually required if SUSY is to solve the naturalness problem associated 
with the Higgs boson mass. They have been the subject of several dedicated searches, listed as ‘third- 
generation’ in Table 1 . Within the 19-parameter pMSSM the masses of these third-generation squarks are 
controlled by different parameters than for the first two generations, allowing the t\,i and masses to 
differ subsfanfially bofh from one anofher and from fhe squarks of fhe firsf fwo generafions. 

Figure 6(a) shows fhe fractional exclusion of model poinfs as projecfed onfo fhe plane of fhe mass of 
fhe lighfer of fhe fwo fop squarks and fhaf of fhe LSR As discussed in Secfion 3.3, fhere are relatively 
few model poinfs af low fop squark mass. This is because a large m{t) (and/or a large frilinear coupling 
At) is required fo obfain fhe large quanfum correcfions needed fo obfain fhe observed Higgs boson mass. 
Despile fhere being relatively few poinfs in fhe inifial sample wifh small fop squark masses, one observes 
fhaf when m{t) is below abouf 600 GeV, mosf poinfs are excluded by ATLAS analyses. 

The sensifivify fo direcf producfion of fop squarks can be seen by considering fhe ATLAS exclusion using 
only fhose analyses from Table 1 fhaf largel direcf producfion of Ihird-generalion squarks. The resulfs, 
presenfed in Figure 6(b), show fhaf when only fhese fhird-generalion ATLAS searches are considered, 
good sensifivify confinues fo be observed for a lighfesf lop squark mass up fo abouf 700 GeV. 

A reasonable correspondence is found belween fhe sensifivify fo fhe pMSSM poinfs and fhose of fhe 
simplified-model decay considered in Ref. [57], in which fhe fop squark was assumed fo decay wifh 
cerfainly fo f -i- (including off-shell lop decays). For larger fop squark masses, in fhe range 600 GeV 
fo 700 GeV, Ihis example simplified-model line exfends info a region where fhe observed fraction of 
pMSSM poinfs fhaf are excluded is less fhan 100%. This can be understood, since in Ref. [57] if was 
shown fhaf fhe sensifivify of fhe t search analyses depends on fhe branching ralio of fhe fop squark fo fhe 
LSR When fhe decay proceeds via fhe Iwo-slep process t\ b + x\ followed by wifh 

fhe assumption fhaf mix"^) = m(x^), fhe 95% CL exclusion limif did nol extend beyond m(i\) = 540 GeV. 
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Figure 6: Fraction of pMSSM points excluded in the plane in various cases. The top two plots show the full 
model set, with impact of all searches on the left and only the third-generation searches on the right. The bottom 
row of plots separates the models according to whether the fi is either mostly (c) left handed or (d) right handed. 
There are relatively few pMSSM points at low fi mass for the reasons described in the text. The simplified-model 
limit overlaid [56, 57] is set assuming directly produced top squark pairs, with each decaying to a top quark and 
neutralino, fi ^ The colour scale is as described in Figure 3. 
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(a) All searches (b) Third-generation searches 

Figure 7; The fraction of pMSSM points excluded by (a) all the listed ATLAS Run 1 searches and (b) just the third- 
generation searches in the b-x\ mass plane. The white line shows a simplihed-model limit [69] made assuming 
directly produced bottom squark pairs, with bottom squarks decaying to a bottom quark and a neutralino, b b+x^y 
The colour scale is as described in Figure 3. 


Since a range of such branching ratios is found in the pMSSM model sample, it is to be expected that 
there is partial sensitivity in this intermediate region. 

The lower two plots in Figure 6 again show that the sensitivity of the ATLAS analyses can depend in a 
non-trivial way on sparticle properties that are not captured by simplified models. The dependence of the 
sensitivity on the left versus right chirality of the top squark is caused in part by different branching ratios 
for decays both to and from that squark. The branching ratios of decays involving wino-like gauginos 
are affected by the SU(2) coupling of the top squark, resulting in significant differences in fhe sensifivify 
depending on whefher fhe ti is dominanfly leff- or righf-handed. Furfhermore, if fhe t\ is mosfly leff- 
handed fhen fhe fop squark’s SU(2) parfner, fhe b, would have a similar mass, allowing analyses fargefing 
sboffom producfion fo become relevanf in consfraining fhe models. 

Figure 7(a) shows fhe sensifivify of ATLAS analyses as projected onfo fhe mass plane of fhe lighfesf 
boffom squark versus fhe LSP. The sensifivify is generally well capfured by a simplified model confaining 
only an addifional boffom squark and fhe LSP [69] . When searches for sparficles ofher fhan fhose of fhe 
fhird generafion are omiffed (Figure 7(b)) fhe similarify of fhe pMSSM wifh fhe simplified model becomes 
still clearer. The sensifivify of fhe ATLAS searches was found fo be similar regardless of whefher fhe 
boffom squark was dominanfly leff- or righf-handed. 


5.1.3. Electroweak sparticles and sleptons 

The sensifivify fo selecfrons, smuons, and fheir sneufrino counferparfs (here denoted collectively by {) is 
displayed in Figure 8. If can be seen fhaf fhe ATLAS searches have good sensifivify fo slepfon masses up fo 
abouf 200 GeV, parficularly when fhe LSP mass is lighfer fhan abouf 75 GeV, where only bino-like LSPs 
survived fhe preselecfion of Table 3. Nevertheless, fhe reduced sensifivify in fhe near-degenerafe region 
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Figure 8; The fraction of pMSSM points excluded by just the electroweak ATLAS searches listed in Table 1, 
projected onto the mass plane, where (is the lightest slepton of the first two generations. The white line reflects 
one of the simplified-model limits (Figure 8(c) of Ref [54]) made for direct slepton pair production assuming that 
left- and right-handed selectrons (or smuons) are mass-degenerate and that each decays via The 

colour scale is as described in Figure 3. 


where the slepton decay produces soft leptons and small means that existing ATLAS searches cannot 
place any lower bound on the slepton mass. The region of sensitivity in the pMSSM projected plane is 
found to have a degree of correspondence with one of the simplified models of Ref. [54]. Generally, 
reduced sensitivity is found in the pMSSM when compared to the more-constrained simplified model. 
This can be understood by recognising that this particular model presupposes that the left- and right- 
handed selectrons and smuons are all mass degenerate, and that each has a 100% branching ratio to a 
lepton and a LSR Breaking these assumptions reduces the number of signal events, and hence allows 
more models to evade detection. 

When the assumption of degenerate left- and right-handed states is dropped from the simplified model, 
the resulting limits are similar to those of the pMSSM. This can be seen in Figure 9, showing the pMSSM 
space projected separately onto the mass of the left-handed or right-handed slepton. The fraction of model 
points excluded is compared to simplified models in which either only left- or right-handed sleptons (e 
and ju) are produced. ATLAS searches have more sensitivity to the production of left-handed sleptons 
since right-handed states lack SU(2) couplings and so have a smaller production cross-section. 

Figure 10 shows the fraction of model points excluded by just the electroweak ATLAS searches listed in 
Table 1, projected onto the plane of LSP and lightest stau mass. It can be seen that the Run 1 sensitivity 
to staus is limited, with large fractions of model points surviving even at the lowest stau masses. This is 
largely because it is difficult to trigger on events resulting from direct stau production, and backgrounds 
to stau searches are much larger than for the equivalent search for sleptons of the first two generations. 
No definitive lower bound can be placed on the stau mass by ATLAS after Run 1. 

Figure 11(a) shows the fraction of models excluded by only the electroweak ATLAS searches (see 
Table 1), this time projected onto the plane of the masses of the lightest two neutralinos. Two prom¬ 
inent features are visible. For x^^ masses lighter than about 200 GeV, a large fraction of models are 
excluded, particularly as becomes large. The dominant exclusion mechanism for large mix^) is 
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Figure 9: Impact of electroweak searches (as listed in Table 1) in (a) the 4-Ti (b) 4 -Lj projections. It 
should be noted that in the 19-parameter pMSSM, the first- and second-generation sleptons of each handedness are 
required to be degenerate. The simplified-model limit in the 4 (4) case is set assuming directly pair-produced left 
(right) handed selectrons/smuons, decaying to an electron/muon and neutralino. The simplified-model limits are 
from Figures 8(a) and 8(b) of Ref. [54]. The colour scale is as described in Figure 3. 
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Figure 10; Impact of electroweak searches (as listed in Table 1) on the f-x^ plane. The colour scale is as described 
in Figure 3. 
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(a) Neutralinos (b) Chargino-neutralino 

Figure 11; Impact of electroweak searches (as listed in Table 1) (a) on the^ 2 “/fi plane and (b) on the^^-jp® plane. 
The 95% CL observed exclusion limit from Ref. [54] is for a simplibed model that assumes pure-wino xf + X 2 
production, followed by the decays ^ W*x^Z*Xy The colour scale is as described in Figure 3. 


due to the Disappearing Traek analysis and is strongest when the LSP is wino-dominated. The gaugino 
mass difference Am^ = m(xj) - m(x^) is typically less than a few hundred MeV for winos and of order 
a few GeV for Higgsinos. As the x^ mass decreases, approaching that of the x^, there is more neutralino 
mixing, leading to a larger Am^, and a shorter^* lifetime, hence the Disappearing Track analysis loses 
sensitivity. The Figure 11(a) row in which m(x^) ~ 50 GeV has lower sensitivity for the Disappearing 
Track analysis. This region is dominated by models for which the relic density is controlled by the Z and 
h boson funnels, so has bino-like LSPs with a Higgsino admixture. Such models do not typically feature 
long-lived charginos. 

For mix^) ^ 400 GeV and mix'^) < 200 GeV, direct production of (and/or jf*) states provides sens¬ 
itivity via the 2-leptons, 3-leptons and 4-leptons analyses. The sensitive region for these multi-lepton 
analyses is similar to that shown from the simplified model of Ref. [54]. Neverfheless fhere remain many 
viable pMSSM poinfs wifhin fhe region excluded in fhe simplified-model scenario. For example, many 
poinfs in fhe Z and h boson funnel regions (m(x^) ~ 50 GeV) have liffle sensifivify in fhe mulfi-lepfon 
analyses as fhe x^ is predominanfly Higgsino-like, leading fo a lower production cross-secfion. 

The equivalenf plof for fhe projection onfo fhe plane of fhe lighfesf chargino and fhe LSP is shown in 
Figure 11(b), again showing fhe fraction excluded by fhe elecfroweak ATLAS searches. In Ibis figure 
fhe Disappearing Track analysis has sensifivify fo models wifh wino-like LSPs which lie close fo fhe 
leading diagonal where m(x^) is only a liffle larger fhan m(x^). Models wifh Higgsino-like LSPs also 
lie close fo fhaf diagonal, buf have larger mass spliffings and so liffle sensifivify from fhe Disappearing 
Track analysis. Away from fhaf diagonal only bino-dominafed LSPs are found. Here fhe besf sensifivify 
is from fhe mulfi-lepfon elecfroweak search analyses (2-lepfons, 3-leplons and 4-lepfons), parficularly 
for mix^) ^ 400 GeV and m(/y°) ^ 200 GeV. The region wifh sensifivify fo fhe mulfi-lepfon searches 
again shows some similarify wifh fhe simplified-model limif from Ref. [54], buf again no region is fofally 
excluded. 
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Figure 12: The fraction of model points excluded by ATLAS Run 1 searches, including those for heavy Higgs 
bosons, projected into the m^-tan^S plane. The white line overlaid is the observed limit for the scenario from 
the HjA TT search of Ref. [75] where the region above the line is excluded. 


5.1.4. Long-lived squarks, gluinos and sleptons 

As described in Section 4.2, model points with long-lived squarks, gluinos and sleptons are treated sep¬ 
arately as the sensitivity to these models of the regular searches for prompt decays is difficult to assess. 
Instead only the Long-lived particle searches [73, 74] are considered. These analyses have very good 
sensitivity to these model points as can be seen from Table 6. Particularly good sensitivity is found in 
cases where the bottom squarks are long lived. The sensitivity to gluinos is quite poor as these are usually 
too short-lived to be picked up by the long-lived particle searches. 


Long-lived 

Bino LSP 

Wino LSP 

Higgsino LSP 

Particle 

Models 

Excluded 

Models 

Excluded 

Models 

Excluded 

g 

899 (5.2%) 

5.1% 

58 (3.4%) 

3.4% 

9 (0.0%) 

0.0% 

bi 

1252 (99.6%) 

76.4% 

51 (100.0%) 

78.4% 

67 (100.0%) 

80.6% 

h 

345 (56.8%) 

36.5% 

6 (100.0%) 

66.7% 

17 (82.4%) 

47.1% 

Tl 

406 (100.0%) 

31 A% 

2 (100.0%) 

0.0% 

41 (100.0%) 

14.6% 


Table 6: Number of model points with long-lived particles and their exclusion fraction. The percentages in paren¬ 
thesis are the fractions of these model points where the long-lived particle has a lifetime long enough to traverse 
the entire detector. 


5.1.5. Heavy neutral Higgs bosons 

Figure 12 shows the fraction of model points excluded by the ATLAS Run 1 searches (which include 
searches for heavy Higgs bosons) when projected into the m^-tan/? plane. The white line overlaid is the 
observed limit from the ATLAS search [75] for heavy neutral Higgs boson(s), H or A, decaying to tt, as 
interpreted in the scenario. A close correspondence can be observed between the limit obtained in 
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that scenario and the region of the pMSSM space excluded by the Run 1 searches in the general pMSSM- 
19 space. The scenario is seen to give a slightly conservative limit compared to the pMSSM space 
as it has light electroweakinos leading to a lower branching fraction for A/// ^ rr than for most pMSSM 
points with similar m^AanyS values. 


5.1.6. Complementarity of searches 


H/A— 

Long-lived sparticles 
Ih, electroweak 
Disappearing Track 
4-Teptons 
3-l^tons, electroweak 
2-taus, electroweak 
2-leptons, elecj^gweak 
tb-i-E^ ,stop 
Stop with Z boson 
Monoiet, stop 
2b-jets 
2-leptons, stop 
1-lepton, stop 
0-lepton, stop 
t(t/I) + jets 
SS/3-leptons -i- jets 
Monojet 
1-lepton -r jets 
0/1-lepton + 3 b-jets 
0-lepton -I- 7-10 ets 
0-lepton + 2-6 jets 


ATLAS 


36 8 
42 9 
0 0 
55 17 
34 14 
53 25 

75J34 

47 13 


11 8 10 
15 8 12 
0 0 0 

12 14 24 
11 35 13 
15 56 19 
26 58 36 

13 17 19 


61 13 14 19 42 
51 42 21 31 34 
30 62 33 43 20 
10 6 16 14 18 



74 

83 

75^ 

39 

35 

40 


4 

48 



44 

\K 

86 


33 

29 



25 

20 


4 

6 

43 

14 

1 

23 


58 

75 

46 

59 

20 

21 


55 

2 


|58 

29 

25 

62 

27 

65 

27 

30 

43 1 


|10 



83 

M 

29 

30 

34 



5 

27 

17 1 

57 

68' 

52 

37 

17 

21 


[57 

2 

32 

91 

'73 

64 

■Ml 

32 

51 


[^ 

71 

7 

22 

79 

72 

68 

88 1 

29 


63 

69 

77 

20 

25 


40 

23 

21 


7 

9 

35 

18 

2 

18 


61 

68l 


26 

26 

32 

61 

67 

6 

18 

SB 

57 


[AT 

26 

18 

21 

73 

58 

2 

26 



[64 

63 

50 

22 

28 

69 

56 

2 

23 


24 

24 

23 

30 

6 

8 

25 

18 


8 



0 14 
25 4 65 

7 19 22 

8 7 58 
15 4 44 

7 12 56 
5 16 61 
35 2 16 

8 5 39 

9 4 39 
23 5 39 

10 1 11 


CD 


"2 O O 

■ 4 —' • 4 —' - 

.S!2, c/3 CD 


Q. Q. ■ 
CD CD 


— CD 


CM 


31 0 
28 0 


13 0 
25 0 


2 1 
10 10 


28 0 
20 0 
16 0 
25 0 
19 0 


O CO CO 
■K CD CD 

w O O 
^ 1-^ ^ 
LU O O 
+ _o _o 
^ CD CD 


CD 


CO 


CO 


Figure 13: Complementarity of ATLAS searches: the figure shows the percentage of model points excluded by the 
analysis on the //-axis that were also excluded by the analysis on the x-axis. References for the individual analyses 
can be found in Table 1 . As none of the models considered are excluded by the ih analysis, it has no overlap with 
any other search. 

The ATLAS searches are designed to be sensitive to different final states, or different kinematic regions. 
Nevertheless it is often the case that a model which produces an excess in one search can have a significant 
signal expectation in others, since related cascade decays can lead to complementary final states. The 
degree of complementarity is explored in Figure 13, which shows the fraction of model points excluded 
(at 95% CL) by one analysis that are also excluded by another. For example, one can observe that of the 
model points that were excluded by the 1-lepton -i- jets -i- F"™®® analysis (fourth row from bottom), 93% 
were also excluded by the 0-lepton -i- 2-6 jets -i- F"™®® analysis (first column). 

Figure 13 demonstrates that there is a good complementarity - relatively small overlap in sensitivity - 
between searches for strongly interacting particles, which are characterised by final states with jets, and 
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Stop with Z boson 

0.4% 

1.0% 

0.4% 

0.5% 

stop 

4.2% 

1.9% 

3.1% 

5.0% 

bh, electroweak 

0 

0 

0 

0 

2-leptons, electroweak 

1.3% 

2.2% 

0.7% 

1.6% 

2-t, electroweak 

0.2% 

0.3% 

0.2% 

0.2% 

3-leptons, electroweak 

0.8% 

3.8% 

1.1% 

0.6% 

4-leptons 

0.5% 

1.1% 

0.6% 

0.5% 

Disappearing Track 

11.4% 

0.4% 

29.9% 

0.1% 

Long-lived particle 

0.1% 

0.1% 

0.0% 

0.1% 

HI A T+T- 

1.8% 

2.2% 

0.9% 

2.4% 

Total 

40.9% 

40.2% 

45.4% 

38.1% 


Table 7; Percentage of model points excluded by the individual analyses. It should be noted that the fraction of 
model points that can be excluded will depend on the model employed and range of input masses initially generated. 
The reader is reminded (Table 2) that the sparticle mass terms in this paper extend to 4TeV. References for the 
individual analyses can be found in Table 1. 
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Figure 14: The density of pMSSM points projected onto the plane of dark matter relic density versus LSP mass, 
before and after the constraints from the search analyses. The colours labelling the different LSP types, as defined 
in Table 4. 


searches for electroweak production. Further study shows that, for the sampling of pMSSM points made 
in this paper, the analyses with the largest regions of unique sensitivity are the 0-lepton -i- 2-6 jets -i- F"™®® 
analysis [58], and the Disappearing Track analysis [72]. Nevertheless some care is required in interpreting 
these results. The degree of apparent overlap is subjective, in that it depends, in some cases sensitively, on 
the metric used when sampling the pMSSM space. Even in cases where the apparent overlap appears to 
be large, for example between the 0-lepton -i- 2-6 jets -i- E™'*® and 0-lepton -i- 7-10 jets -i- E™®® analyses, 
both searches are found to have regions of pMSSM space in which they provide unique sensitivity. The 
Disappearing Track analysis is mostly sensitive to model points with a wino-like LSP, so an alternative 
prior (or weighting by LSP type) of the sample model points would directly affect the apparent relative 
sensitivity of this analysis. 

The overall fraction of model points within the pMSSM space excluded by each analysis for each of 
the LSP types is shown in Table 7. Only the €h analysis is unable to constrain the pMSSM set with 
the luminosity available. The lack of sensitivity for that analysis is not unexpected since for simplified 
models it excludes only points with very light LSPs [70]. It should again be noted that the absolute 
values of the fractions of model points excluded is strongly affected by the prior sampling, in particular 
by the upper mass bounds used for the scan in selecting the pMSSM input parameters (see Table 2). 
The relative fractions of model points excluded by each analysis are a little more informative, but again 
care is necessary in their interpretation since they too are sensitive to changes to the assumptions or 
constraints applied to the initial model set. Nevertheless, the high sensitivity of the 0-lepton -i- 2-6 jets -i- 
analysis for all LSP types, and the Disappearing Track analysis for models with a wino-like LSP is 
unambiguous. 
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5.2. Impact of ATLAS searches on dark matter 

The nature of the LSP has a strong influence on the expected dark matter relic density. For the pMSSM 
points, the initial mass spectrum of the LSP - before applying any ATLAS SUSY search constraints - is 
sculpted by the requirement that the dark matter relic density, as predicted from thermal production and 
annihilation calculations in the early universe, should not be larger than that observed. As discussed in 
Section 3.2.2, the LSP is not required to account for all of the dark matter density, since other particles 
may contribute. 

The effect of the relic density preselection can be seen in Figure 14(a), which shows the density of 
pMSSM points in the plane of the^j relic density (Q.(x^)h^) versus the mass of the^®. The model points 
with a bino-like LSP are shown in red, while those with wino-like and Higgsino-like LSPs are in blue and 
green respectively. The features at low LSP mass are due to the effective annihilation of LSPs through 
5-channel Z or Higgs bosons - the so-called Z-funnel and Higgs-funnel regions. There are few wino- 
or Higgsino-dominated LSPs at low mass since in such cases the is expected to be accompanied by 
an almost degenerate chargino, which would have been observed at LEP [102]. Most of the models with 
wino- and Higgsino-dominated LSPs lie on bands which are almost straight lines on the logarithmic plot, 
since the thermally averaged annihilation cross-section is expected to be proportional to the inverse square 
of the LSP mass, resulting in almost exact proportionality between and 

Figure 14(b) shows the distribution of model points in the same plane after applying the constraints from 
the ATLAS SUSY searches. Considering first the bino-dominated LSP model points, one observes that 
model points are excluded by the ATLAS searches across a wide range of LSP masses (and expected relic 
densities). It is found that the ATLAS searches have some sensitivity up to m(x^) < 800 GeV, while at 
low bino mass about two-thirds of the LSP model points in the Z- and Higgs-funnel regions are excluded. 
For wino-dominated LSPs, the overall range of sensitivity in m(x^^) again extends up to about 800GeV. 
Model points with m(y^j) ^ 220 GeV are particularly depleted by the ATLAS searches with about 80% 
of model points with wino-like LSPs in this mass range being excluded, mostly by the sensitivity of the 
Disappearing Track analysis to the charged wino. For Higgsino-dominated LSPs, the ATLAS sensitivity 
is smaller as thex^-x^ mass splittings are mostly too large to have an observable^* lifetime. 

Focusing on model points with bino-dominated LSPs, Figure 15 shows the model point density before and 
after ATLAS Run 1 searches. In this figure the colour code now corresponds to the dominant annihilation 
mechanism of the dark matter - for example an orange point has at least one squark of the first two 
generations, or the gluino, close in mass to the LSP, allowing annihilation mechanisms such as q +x^ —> 
q+y to proceed effectively. The vertical cut-offs for each coannihilafor correspond fo fhe preselection from 
previous experimenfal constrainfs - for example charged sparficles lighfer fhan 100 GeV are forbidden. 
The differenf minima seen for fhe relic dark matter energy densify for each coannihilafor correspond fo 
fhe differenf coupling sfrengfhs. 

The ATLAS searches are seen fo be particularly sensifive fo model poinfs wifh lighf-flavour squark or 
gluino coannihilafors. One would expecf such model poinfs, which have af leasf one coloured parficle 
near-degenerale wifh the LSP, to produce events with missing transverse momentum, and soft jets from 
the coloured sparticle decays. Such model points are difficult to observe in the inclusive analyses, but 
when produced in association with ISR can be observed through the monojet analyses [64, 67]. 

The limits from experiments searching for direct detection of dark matter were taken into account when 
generating the pMSSM points. Figure 16 shows the corresponding predictions together with the limits 
set by the measurements of the spin-independent (SI) and the proton spin-dependent (SD) cross-sections 
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Figure 15; The density of pMSSM points on the plane of relic density versus LSP mass, before and after the 
ATLAS constraints, for the bino-like LSP models. The model points are colour coded by annihilation mechanism 
as described in the text. ‘Gaugino’ refers to electroweak gaugino coannihilation, while ‘Light Flavour’ refers to 
coannihilation either against squarks of the first two generations or against gluinos. 


by the LUX [98] and COUPP [99] experiments respectively. The direct detection limits include the 
additional uncertainties described in Section 3.2.2. In each case the interaction cross-section is scaled 
by a factor of Rq = Q.{x^^)h^IQ.Y{aackh^, since in this paper it is assumed that the LSP may be only one 
of a range of possible contributors^ to the dark matter abundance, whereas direct detection experiments 
interpret their results in a framework in which the LSP completely saturates the relic density. 

Figure 16(a) shows the model points projected onto the plane of the spin-independent interaction cross- 
section versus the LSP mass. The direct detection limit from the LUX experiment is shown, in which, 
as already mentioned, the relic density of the colliding dark matter was assumed to saturate the value 
measured by the Planck Collaboration. The only model points with LSP mass less than 100 GeV are those 
with bino-like LSPs lying in the Higgs- or Z-funnel regions. Around 15% of the bino-like model points 
that pass the preselection (Table 3) would have been excluded by the nominal LUX limit had the additional 
uncertainty scaling factor of four, as mentioned in Section 3.2.2, not been applied. Figure 16(b) shows 
the same space after removing model points excluded (at 95% CL) by the ATLAS Run 1 searches. 

Figure 16(c) shows the spin-dependent cross-section, again as a function of the neutralino mass. The plot 
also shows the corresponding direct detection exclusion limit, this time from the COUPP experiment. The 
same plot is shown after ATLAS Run 1 search constraints in Figure 16(d). Again one can see that sens¬ 
itivity from LHC searches can stretch to regions with cross-sections several orders of magnitude below 
those of the current best direct detection experiments. Similar conclusions are found for spin-dependent 
neutron cross-sections, in this case the best direct detection limit being from Xenon-100 [100]. 


Figure 16(b) shows a roughly rectangular region excluded by ATLAS searches, bounded approximately 

-0\ ^ d- ^SI 

N-x\ 


by m(jLj) 220 GeV and x 10 A similar region is excluded in the SD case (Fig¬ 


ure 16(d)), this time for R^i x .q ^ 10 for the same LSP mass range. These regions are 
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® These other possible contrihutors are not descrihed hy the pMSSM points considered in this study. 
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Figure 16: Left, the distribution of model points generated; right, the distribution of model points not excluded by 
ATLAS Run 1 searches, as projected onto the scaled spin-independent (SI) interaction cross-section of nucleons 
with the neutralino versus the neutralino mass. The cross-sections are scaled by a factor of Rq = f2(Jj)/i^/f2pianck/i^- 
The calculated spin-independent interaction cross-sections are a weighted average of the contributions from proton 
and neutron scattering, corresponding to the Xenon atom (the target nucleus of the LUX experiment) and normal¬ 
ised to one nucleon. The 90% confidence limit [98] from the LUX direct detection experiment is overlaid, in which 
it is assumed that the dark matter comprises only the LSP, with relic density as measured by the Planck Collabora¬ 
tion [97]. For the spin-dependent cross-sections, the calculated proton cross-section is shown. It is compared to the 
direct detection limit from the COUPP experiment [99]. 
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Figure 17; Distributions of pMSSM points before and after the constraints from the ATLAS direct searches for 
observables related to the decays of the observed Higgs boson, (a) Higgs-h-quark coupling ratio Kb, for all LSP 
types, (b) Branching ratio of the Higgs boson to the LSP, only for those pMSSM points with bino-like LSPs and 
< 65 GeV. The dashed lines show the observed and expected upper bounds of the 95% conhdence limits on 
those parameters. The subplots show the fraction of model points excluded by the ATLAS direct searches. 


dominated by models with an LSP that is almost pure wino. The winos, having only very small Higgsino 
admixtures and thus little interaction with the Higgs boson, have small direct detection cross-sections. 
Such model points then have small mass splittings, and are tightly constrained by the Disappear¬ 

ing Track analysis, since the lightest chargino (the charged wino) is sufficiently long-lived to be detectable 
in that analysis. As the degree of Higgsino mixing increases, the direct detection cross-sections increase, 
and so does the LSP-NLSP mass splitting. This results in a more rapid NLSP decay, and lack of sensit¬ 
ivity by ATLAS. 

As one would expect, the direct detection and LHC search techniques are found to be very complementary. 
The ATLAS searches have sensitivity for m(x^) up to about 800 GeV, and are particularly effective in 
constraining low-mass LSPs. ATLAS shows sensitivity to many SUSY models with direct detection 
cross-sections several orders of magnitude below current direct detection limits. On the other hand, direct 
detection experiments rule out large x cr up to higher LSP masses than the LHC can access. 


5.3. Effect of ATLAS Higgs boson coupling measurements 

Direct searches for SUSY particles are complemented by indirect searches for the effects of new particles 
and fields. In SUSY models the presence of the second Higgs doublet generally results in light Higgs 
boson couplings that are modified with respect to their SM values. For example, enhancements to the 
couplings to down-type quarks and charged leptons are predicted at large tan/3. In addition, loop con¬ 
tributions from virtual SUSY particles can produce an observable effect on lower-energy observables, 
including Higgs boson decay rates. 
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Combined fits by ATLAS to the measured Higgs boson production and decay rates provide a sensitive 
probe for such additional Higgs fields and for SUSY particles contributing either through loop effects or 
through invisible decays of the Higgs boson to unobserved neutralinos. The relevant ATLAS analyses, 
statistical methods, and systematic uncertainties are described in Refs. [143, 144]. These measurements 
are based on up to 4.8 fb“' of pp collision data ai = 1 TeV and up to 20.3 fb“^ at = 8 TeV. They 
are used to constrain coupling scale factors 


Ki = 



BR(/j —> i + i) 
BR(/j ^ i + i)sM ’ 


where T/, represents the total width of the observed Higgs boson, BR indicates its branching ratio, and 
a SM subscript indicates the Standard Model prediction. The measurements of these scale factors were 
derived using the observed Higgs boson’s visible decay channels as reported in Ref. [144]. The upper limit 
on the invisible branching ratio, taking account of both visible and invisible Higgs boson decay channels 
is reported in Ref. [145].^ Most of the coupling parameters are not yet sufficiently well measured to 
have an effect on the SUSY parameter space. The exceptions - those two parameters which are already 
sufficiently well measured to constrain the pMSSM - are discussed below. 

The parameter Kh corresponds to the coupling of the Higgs field to t^-quarks. Measuring a value of Kh 
differing significantly from unity would be a clear indication of physics beyond the SM, and an important 
constraint on the allowed pMSSM parameter space. For the current ATLAS fit, the expected 95% CL 
range is \Kb\ < 1.8. This expected region contains almost all the pMSSM points not excluded by the direct 
searches, so the expected constraint from Kb is negligible. The current ATLAS observed fitted value for 
Kb is 0.62 ± 0.28 at 68% CL [145]. This value is lower than, but consistent with, the SM value. The 
corresponding observed 95% CL spans the range \Kb\ < 1.2, and is more restrictive than the expected 
range. As a result the observed fit disfavours those pMSSM points (see Figure 17(a)) with Kb > 1.2, 
which typically are those with large tan/3. Those models that are disfavoured at the 95% CL represent a 
weighted fraction of 3.1% of all the pMSSM points. This fraction is rather similar before and after the 
ATLAS Run 1 direct searches are considered. 


The observed Higgs boson can decay to a pair of LSPs if those LSPs are sufficiently light. This can 
lead to an enhanced branching ratio of the h to the invisible final state. The expected branching ratio for 
decay of a light Higgs boson via ZZ* to neutrinos is (9(10~^), much smaller than the present experimental 
uncertainties, and so is negligible in this context. No models with a sufficiently light wino- or Higgsino- 
like LSP satished the initial constraints of Section 3.2, so only models with a bino-dominated LSP need 
be considered in this context. All such LSPs have some wino and Higgsino admixture, allowing them to 
couple to the Higgs boson. Figure 17(b) shows the calculated distribution of the lightest Higgs boson’s 
invisible branching ratio for those pMSSM points with bino-like LSPs and having mf/fj) < 65 GeV. To 
account for hnite-width effects, this selection allows for models with LSP masses a little above m{h)l2. 
These light-LSP models represent 4.5% of the model points with bino-like LSPs that survive the ATLAS 
searches in Table 1. The observed ATLAS bound on the h invisible branching ratio is BR(/j ^ +x^) < 

0.23 at the 95% CL [145]. This observed value is close to the expected upper limit of 0.24. Considering 
only those pMSSM points with a bino-like LSP lighter than 65 GeV, and which were not excluded by 
the direct searches, the fraction disfavoured at the 95% CL by the invisible branching ratio ht is 6.4%. 
The disfavoured points correspond to a weighted fraction of 0.0056% of all the pMSSM points (with any 
LSP type) not excluded by the ATLAS direct searches. The corresponding expected fractions are 6.1% 

^ These constraints are based directly on the log-likelihood test statistic, rather than a calculation of CL, hence these particular 
limits are not protected against disfavouring models where there is little sensitivity. 
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(of those surviving models having bino-like LSPs with < 65 GeV) and 0.0053% (weighted by LSP 
type) of all surviving models. 

There is no overlap between the models disfavoured by the Ky measurement and those disfavoured by 
the Higgs boson invisible branching ratio. This demonstrates that these indirect searches for new physics 
complement the direct searches in different ways by constraining different parts of the pMSSM space. 


5.4. Impact of ATLAS searches on precision observables 

Loop contributions from SUSY particles can also affect lower-energy observables. The measured val¬ 
ues of several precision observables were taken as prior constraints on the pMSSM space sampled, as 
explained in Section 3.2.1. 

The effect of the ATLAS Run 1 searches on the distribution of pMSSM points as projected onto the 
expected values of these precision observables is shown in Figure 18. In each case the most noticeable 
feature is that the ATLAS direct searches remove pMSSM points rather uniformly across the space of each 
precision variable, demonstrating the complementarity of the searches. The direct searches from Run 1 
have placed tight constraints on the strongly interacting sector of the pMSSM, whereas the precision 
variables depend on other sparticles, for example the smuon and electroweakino masses, which are less 
tightly constrained by direct searches so far. 

Figure 18 supports the statement in Section 3.2.3 that the exclusion power of the ATLAS searches does 
not depend significantly on the mass range used for the light (SM-like) Higgs boson, the fraction of 
excluded model points being remarkably flat. The other three distributions: BR(B^ ^ yu/i), BR(f7 ^ ^y), 
and the A(^ - 2)^ also show good exclusion by ATLAS for all values of the observables. It is noticeable 
that the majority of the pMSSM points have only small SUSY contributions to these observables, so 
those model points would not be expected to be discoverable by using those measurements in the near 
future. Nevertheless there exists, for each precision measurement, a small number of model points with 
SUSY contributions large enough to be discovered by each. One observes that ATLAS also tends to have 
slightly larger sensitivity to these ‘tail’ model points. This increased sensitivity is not surprising since 
these model points are generally those with some lighter SUSY particles contributing to the precision 
measurements via loop diagrams. This correlation becomes more clear in Figure 19, which shows the 
number of generated model points in the plane of the K{g - 2)^ versus the mass of the lighter left- or 
right-handed smuon (min(m(/iL), w^Our)))- The experimentally measured value [89], represented by the 
hatched band, only overlaps a region where m{fi) < 1 TeV. As noted previously in Section 3.2.1, the 
allowed range used for A(g - 2)^ in this paper is the union of the 3cr intervals around the SM value and 
the experimental measurement. The experimentally measured value, if confirmed, would be a powerful 
consfrainf on fhe space. 


5.5. Fine-tuning 

The nafuralness of pMSSM poinfs is nof considered in fheir generation, buf if is inferesfing fo see jusf 
how fine-funed fhe model poinfs are affer fhe ATLAS exclusion. The fine-funing paramefer adopfed here 
is fhe one defined by Barbieri and Giudice [146]. Figure 20 illusfrales fhe prediction of fhe fine-funing 
before and affer ATLAS exclusion, for all LSP fypes (a), and for each LSP-fype separafely (b)-(d). The 
shapes of fhe disfribufions are determined by fhe ranges of fhe pMSSM parameters considered. The 
dominanf confribufions fo fhe fine-funing come from /r and At. If can also be observed fhaf fhe ATLAS 
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Figure 18; Distributions of model points before and after applying ATLAS searches for various precision ob¬ 
servables. Figure (a) shows the mass of the light (SM-like) Higgs boson, m(h), (b) shows the branching ratio of 
Bs ^ fj. + fi, (c) is the branching ratio of Z? —> s + j and (d) is the difference between the predicted value of gr - 2 
and the SM value, A{g - 2)^. The subplots show the fraction of model points excluded by ATLAS as a function of 
the observables. 
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Figure 19: Distribution of model points in the plane of the mass of the lightest left- or right-handed smuon versus 
A(gf - 2)^. The experimental measurement is overlaid as the hatched band [89]. 


Run 1 searches exclude pMSSM points with a wide range of different fine-tuning values. A decrease 
in sensitivity is observed with increasing fine-tuning for models with a Higgsino-like LSR This is an 
artifact caused by the fact that the Higgsino is the LSP, that its mass is driven by the ji parameter, and 
so by requiring large fine-tuning, one is indirectly requiring that there be a heavy LSP, reducing ATLAS 
sensitivity. 

The lowest fine-tuning parameter value for a surviving model point is 56. A representation of the sparticle 
spectrum of this point can be found in Figure 21(a). As expected for a generic model point with low fine- 
tuning, the top squarks have sub-TeV masses. The model point has a wino-dominated LSP with a mass 
of 107 GeV, which is lighter by 1.6 GeV than its charged partner. This mass difference is sufficient to 
ensure that the chargino lifetime is short enough (< 1 ps) to evade the Disappearing Track analysis. The 
lower part of the mass spectrum includes a Higgsino-dominated {x\, X^, T 2 1 multiplet with masses around 
230 GeV which are mixed with the winos at about the 10% level. Most strongly interacting sparticles, 
particularly the gluino and left-squark doublet, have large masses of around 3 TeV, which is beyond the 
Run 1 search reach. The heavy Higgs bosons are also beyond the Run 1 search reach. 

Similar model points survive with a Higgsino-like LSR for example as shown in Figure 21(b), which has 
a fine-tuning of 57. Model points such as those in both Figure 21(a) and Figure 21(b) would be expected 
to produce observable signals at LHC Run 2 through a variety of channels. In particular in both cases the 
lighter top squark, with mass of about 800 GeV is not far from the current ATLAS Run 1 sensitivity. 

Aside from model points with a wino-dominated or a Higgsino-dominated LSP there also remain model 
points with low fine-tuning and a mixed LSP. For example. Figure 21(c) shows a model point with a 
mixed bino-Higgsino LSP a so-called ‘well-tempered’ neutralino case [105]. This model has a bottom 
squark with mass of about 650 GeV, not far above the Run 1 direct search sensitivity. 

Models with a small number of low-mass squarks of the first two generations also survive. For example. 
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Figure 20: Distribution of fine-tuning (as defined in Ref. [146]), before and after ATLAS exclusion. The subplots 
show the fractions of model points excluded by the ATLAS Run 1 searches. 
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Figure 21: The sparticle spectra for some of the pMSSM-19 model points with the smallest fine-tuning not to have 
been excluded by the ATLAS Run 1 searches. The dashed lines indicate the dominant decay modes. For more 
details see the text. 
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Figure 21(d) shows a spectrum of a model point with a % squark with a mass around 800 GeV, which 
again would be likely to be discoverable at LHC Run 2. 


6. Conclusion 

The ATLAS Collaboration has performed a wide range of direct searches for supersymmetry during 
the first run of the LHC, using pp collisions with centre-of-mass energy up to 8TeV and an integrated 
luminosity of up to 20.3 fb“^ The interpretation of those results within the wider framework of the 
pMSSM gives insights into the breadth of sensitivity of those searches. 

From an initial random sampling of 500 million pMSSM points, generated from the 19-parameter pMSSM, 
a total of 310,327 model points with^° LSP are selected each of which satisfies constraints from previ¬ 
ous collider searches, precision measurements, cold dark matter energy density measurements and direct 
dark matter searches. The models are importance-sampled so that there are roughly equal numbers of 
models with a bino-like, wino-like or Higgsino-like LSR For these model points, more than 30 billion 
signal events are generated, of which more than 600 million events from 44,559 different model points are 
passed through a GEANT4-based fast detector simulation and full reconstruction to accurately determine 
which have been excluded (at 95% CL) by ATLAS searches. 

The impact of the ATLAS Run 1 searches on this space is presented, showing their overall effect in con¬ 
straining such supersymmetric models. The results are particularly clear when considering the fraction of 
model points surviving, after projection into two-dimensional spaces of sparticle masses. The constraints 
on the masses of squarks (including third generation squarks), gluinos, electroweakinos, sleptons and 
heavy neutral Higgs bosons are all presented. 

For the energy and luminosity achieved in LHC Run 1 the ATLAS constraints are most effective for 
strongly interacting sparticles, with weaker constraints on electroweakinos and sleptons. A general con¬ 
gruence is observed between the pMSSM points excluded and the limits determined previously in the 
context of simplified models. Nevertheless, significant differences are observed depending, for example, 
on the number of kinematically accessible squarks, and on their flavour and couplings, particularly for 
direct sleptons production. 

The most constraining ATLAS analyses - for the model points generated - were the 0-lepton -i- 2-6 jets 
-I- analysis, and the Disappearing Track analysis, the latter being especially powerful in the case of 
model points with a light wino-like LSR Good complementarity is observed between different ATLAS 
analyses, with almost all showing regions of unique sensitivity. 

When considering dark matter predictions, the LHC experiments are very complementary to the direct 
detection experiments. The two techniques have different sensitivities across the pMSSM-19 parameter 
space, with the ATLAS searches having more sensitivity at lower LSR mass (with significant sensitivity 
for nr(Yj) < 800 GeV) and the direct detection experiments providing bounds at larger LSR-nucleon 
scattering cross-section. Similarly, the ATLAS experiment provides constraints which are complementary 
to those of precision measurements of BR{Bs pp), BR(b —> sy), and A{g - 2)^. 

Model points with relatively low fine-tuning (of order 50) remain viable after LHC Run 1. Those with 
the lowest fine-tuning have relatively light top squarks, which indicates that one would expect them to be 
accessible by ATLAS searches with the LHC Run 2 dataset. 
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Appendices 


A. Model point calculation 

The calculation of model point properties proceeds as follows. First, SoftSUSY 3.4.Q [147] is used, 
embedded in micrOMEGAs 3.5.5 [77, 78], to calculate the sparticle spectrum. This spectrum is used 
as an input to FeynHiggs 2.10. Q [103, 104], which recalculates the light CP-even Higgs boson mass 
including important 3-loop corrections not present in SoftSUSY. For sparticle decay tables, a modified 
version of the program SUSY-HIT 1.3 [148] is used, updated to use HDECAY 5.11 [149] for the Higgs 
decays. In some cases, where relevant or even dominant decay modes are not included in the SUSY-HIT 
output, MadGraph5_aMC@NL0 2.1.1 [150] is instead used to recalculate those decays as described in 
detail below. MicrOMEGAs 3.5.5 is used to calculate the dark matter relic abundance along with the spin- 
independent and spin-dependent scattering cross-sections in direct detection experiments. MicrOMEGAs is 
also used to calculate the following flavour physics and precision electroweak observables: Ap, /\{g - 2)^, 
BRff? ^ sy) and BR(Bi ^ Finally, analytic formulas are used to calculate BR(R^ ^ T*Vr) and 

the Z boson invisible width, and to check the stability of the vacuum. 


B. Sparticle decay calculation 


The calculation of sparticle decays is performed as follows. First, the publicly available version of 
SUSY-HIT is modified by: (/) Incorporating light quark and lepton masses in the calculation of branching 
ratios and lifetimes for the various sparticles. Two-body decays implement the full mass corrections, 
while three-body decays only include a modified phase-space cut-off. The mass of fhe lighfesf meson 
of fhe appropriafe fype is included in fhe relevanf phase space calculations fo accounf for hadronisafion 
effecfs. (//) Employing full analyfic expressions from Ref. [151] for chargino decays when fhe chargino- 
neufralino mass spliffing is < 1 GeV, in which case a careful freafmenf of hadronisafion is imporfanf and 
can significanfly affecf fhe model phenomenology, {in) Removing QCD correcfions fo decays involving 
fop and boffom squarks due fo fheir fendency fo resulf in negative decay widfhs. {iv) Separating fhe 
Higgs decays fo r sneufrinos from fhose fo election and muon sneufrinos (fhe public version of HDECAY 
calculafes fhem separafely buf fhen averages fhem). 

Several cases have been identified in which SUSY-HIT does nof calculate imporfanf decay channels or 
predicfs a value fhaf differs significanfly from fhe full mafrix elemenf predicfion, which are defailed in fhe 
following. 


B.l. Right-handed sfermion decays 

For righf-handed sfermions, SUSY-HIT only calculafes fwo-body decays of fhe form fx^ (and 

/r ^ fg for righf-handed squarks). If fhe bino componenf of fhe kinemafically accessible neutialinos is 
very small, fhese decays can be highly suppressed. As a resulf, if is common for fhree-body decays of 
fhe form /r fx^, where B is an electioweak gauge boson, fo dominafe over fhe fwo-body decays. In 
cases where fhe spliffing befween /r and fhe ESP is small, fhe fhree-body decay modes can be forbidden 
and a four-body decay (wifh B off-shell) can pofenfially be dominanf, allhough numerically fhis lurns 
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out to be uncommon. Therefore MadGraph5_aMC@NL0 is employed to calculate the decays of any right- 
handed sfermion satisfying the following selection criteria (designed to select all model points for which 
the multi-body decays would be significant while minimising the number of model points selected): 

1) The bino content of the LSP is less than 10%. 

2) Any accessible neutralinos with a bino content above 10% have a mass splitting with /r that is less 

than 100 GeV or below 20% of the ^-LSP mass splitting. 

3) For squarks, if the gluino is accessible it must have a mass splitting with fn that is less than 100 GeV 

or below 10% of the ^-LSP mass splitting. 

For model points with a - LSP mass splitting above 100 GeV, all possible three-body decays are cal¬ 
culated with MadGraph5_aMC@NL0 and added to the existing two-body decays. This can lead to double¬ 
counting when three-body decay diagrams contain intermediate particles which can go on shell. In this 
case, the on-shell component of the three-body decay is already included in the two-body decay width 
(as a tree-level two-body decay in which the SUSY decay product also undergoes a tree-level two-body 
decay). This is resolved by selectively removing redundant decays. First, the sum of the combined three- 
body decay modes is compared with the sum of tree-level two-body decays for which the SUSY decay 
product has tree-level two-body decays (the “on-shell component”). In the case where the on-shell com¬ 
ponent comprised over 80% of the total three-body width, the three-body decays are discarded and only 
the two-body decays kept, since the sequential two-body decays give a more accurate representation of 
the decay kinematics. If the on-shell component is less than 80% of the total three-body width, the off- 
shell componenf of fhe fhree-body decays is deemed significanl and fhe redundanf on-shell componenf 
(fhe selecfed fwo-body decays) is discarded. 

For model poinfs wifh a /r-LSP mass spliffing below 100 GeV, four-body decays could be imporfanf. 
Here use was made of fhe MadWidfh package in MadGraph5_aMC@NL0 [150] by calling fhe compute_widths 
function; fhis function aufomafically defermines whefher fhree-body and four-body decays are required 
and includes an aufomafic mechanism for removing fhe degeneracy. If is found fhaf MadWidfh reproduces 
fhe resulfs obfained using fhe procedure described above for several fesf cases wifh fhree-body decays. 

In fofal fhe fhree-body decays are recalculafed for af leasf 1 sparficle decay in 160482 model poinfs. 

For 86109 model poinfs, mr decays are recalculafed, while for 86791 model poinfs fhe Jr decays are 
recalculafed and for 122127 model poinfs fhe cr decays are recalculafed. The four-body decays are 
calculafed for 11633 model poinfs, alfhough fhis ends up making a significanf confribufion only for a few 
model poinfs. 


B.2. Wino and Higgsino decays to sfermions 

Alfhough SUSY-HIT normally calculafes fhree-body decays of neufralinos and charginos fhaf have mass 
spliffings wifh fhe LSP below fhe Z or IT boson masses, if does nol calculafe fhese decays when ofher 
fwo-body decay modes are available (fo avoid fhe sorf of double-counfing described above). This can 
be a problem when fhe available fwo-body decays are suppressed by a small mixing angle (which can 
occur if fhe decaying gaugino is a neufral wino or Higgsino and fhe fwo-body decay is fo a righf-handed 
sfermion) or by kinemafics. For fhis reason MadGraph5_aMC@NL0 was employed fo calculafe fhe full 
fhree-body decays in cases where fhe decaying gaugino has a bino confenf below 10% and a righf-handed 
sfermion is accessible, or has an accessible decay fo any sfermion where fhe sfermion-gaugino spliffing 
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is below 20% of the gaugino-LSP splitting. This occurs in 716 model points with wino or Higgsino LSP 
and 2080 model points with bino LSP Once again the issue with double-counting between the two-body 
and three-body decays is resolved; in this case any tree-level two-body decays where the product has a 
tree-level two-body decay are removed, since the off-shell three-body decays are expected to be important 
in nearly all of the selected scenarios. 

B.3. Four-body top squark decays 

Although SUSY-HIT includes a calculation of the four-body top squark decays, a significant discrepancy 
is observed between the SUSY-HIT results and the decay width calculated by CalcHEP in a previous 
model set. Therefore MadGraph5_aMC@NL0 is used to recalculate the top squark decays for any model 
point in which SUSY-HIT predicts a non-zero four-body decay rate for the light top squark. In this case, 
there are no concerns about overlap with other top squark decay modes. This procedure is employed for 
15 model points with wino or Higgsino LSP and 5127 model points with bino LSP. 


C. Calculational pathologies 

To ensure that the spectrum is calculable, any model point with the following SOFTSUSY errors is dis¬ 
carded: “No acceptable solution found”, “Non-perturbative”, “No convergence”, “Inaccurate Higgs mass”, 
“Numerical problemThrown”, and “Not achieved desired accuracy”. Additionally the accuracy of the 
spectrum generation is tested by re-generating the spectrum with SuSpect [32] and ensuring that SuSpect 
does not produce any fatal errors or predict any sparticle masses to differ by more than 50% from the 
SoftSUSY prediction. Since FeynHiggs is used for the light CP-even Higgs boson mass, any model 
point for which FeynHiggs indicates an error or fails to write an output file, is discarded. Additionally 
it is required that m{h) be within 5 GeV of the SOFTSUSY prediction, since a small tail of model points is 
observed for which this deviation is extremely large, suggesting that the calculation may have been un¬ 
reliable. FeynHiggs also gives its own estimate of the uncertainty in the Higgs boson mass calculation. 
When this uncertainty is larger than 5 GeV the model points are discarded as being unreliable. Finally, 
the reliability of the SUSY-HIT decays is asserted by discarding model points for which any particle has 
a width larger than 1 TeV or a negative branching ratio. 


C.l. Theoretical Constraints 

In addition to the numerical pathologies described above, model points that produced the following 
SoftSUSY errors are discarded as being theoretically inconsistent: “tachyon”, “MuSqWrongsign”, “m3sq- 
problem”, and “Higgs potential ufb”. Additionally it is checked that the scalar potential does not break 
colour or charge, using Equation (6) of Ref. [152] for Aj and analogous formulas for A^ and Ar. It should 
be noted that Equation (6) of Ref. [152] uses m(tL) and m{tR), while this paper uses the Eagrangian 
parameters m{Q^) and m{u^). 
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D. Importance sampling by LSP type 


As described in Section 3.3, the model points generated are categorised according to LSP type (see 
Table 4), and importance sampled so that approximately equal numbers of model points are available 
for each LSP type. 

There are various ways in which over-production of dark matter can be achieved. A neutralino LSP can 
annihilate through 5-channel exchange of a Z or Higgs boson or through t-channel exchange of a char- 
gino or sfermion. If another sparticle is nearly degenerate with the LSP, that sparticle can also enhance 
the effective LSP annihilation rate through coannihilation. In general, sparticles with the strongest self¬ 
interactions (particularly coloured sparticles) make the most effective coannihilators. Supersymmetric 
models with a wino-like LSP or Higgsino-like LSP (see Table 4) have a chargino that is nearly degen¬ 
erate with the LSP, resulting in a sizeable annihilation rate through chargino exchange; the degenerate 
chargino also serves as a coannihilator. The relic density constraint is therefore satisfied as long as the 
wino (Higgsino) mass is below 2.7 TeV (1 TeV) [153-155]. On the other hand, bino-like LSPs (A^j > 
max(Aj 2 , -i- have no guaranteed annihilation mechanism. One possibility is for the LSP to mix 
with the Higgsino [156], giving couplings to the Z and Higgs bosons and to the charginos. In order for the 
annihilation rate to be large enough, the Higgsino must be somewhat degenerate with the LSP (resulting 
in a large mixing angle) or the LSP must be approximately half the mass of the Z or MSSM Higgs bosons 
(termed a “funnel” region), producing a resonantly enhanced annihilation rate that compensates for the 
mixing-angle-suppressed couplings. Alternatively, a chargino or sfermion can be nearly degenerate with 
the LSP, resulting in a sizeable t-channel annihilation rate and/or significant coannihilation [157]. 

Correcting the undersampling of model points with bino-like LSPs is important for two reasons. First, 
they are the only model points in which the observed relic abundance can be obtained with an LSP 
light enough to be seen at the 8 TeV LHC, making them particularly interesting from the dark matter 
perspective. Second, the absence of an associated chargino can significantly alter the available decay 
modes. As an example, top squark decays to only occur in models with bino-like LSPs, since the 

top squark can simply decay to bx'^ in wino or Higgsino LSP models. 
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